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occur  when  the  primary  radiating  species  of  the  plume  source  and  the 
primary  absorbing  species  of  the  atmosphere  are  the  same  species  (i.e., 
H?0  and  COo).  This  correlation  effect  is  handled  correctly  by  treating 
the  entire  line  of  sight  from  the  sensor,  through  the  intervening  atmosphere 
and  the  hot  source  region  of  the  plume,  as  a single  highly  inhomogeneous 
and  nonisothermal  optical  path.  A detailed  investigation  is  made  of  the 
ability  of  the  standard  Curtis -Godson  and  the  newly  formulated  Lindquist- 
Simmons  approximations  to  handle  correctly  the  high  degree  of  inhomo- 
geneity of  this  optical  path,  and  the  conditions  under  which  the  Curtis - 
Godson  approximation  is  likely  to  fail  are  determined.  For  these 
conditions,  the  Lindquist-Simmons  approximation  must  be  used,  and  pro- 
cedures for  the  efficient  application  of  this  approximation  are  derived. 

A discussion  of  band  model  parameters  used  for  calculations  within  the 
band  model  is  presented.  The  General  Dynamics  high-temperature  band 
model  parameters  and  parameters  derived  from  the  Air  Force  Cambridge 
Research  Laboratories  (AFCRL)  line  data  compilation  for  atmospheric^ 
absorbing  species  are  compared  and  analyzed  to  determine  the  spectral 
and  temperature  regions  for  which  they  are  most  applicable. 

A brief  description  is  given  of  the  computer  code  (ATLES)  that  handles 
the  band  model  computations.  In  addition,  a detailed  analysis  is  made  of 
the  radiative  transfer  properties  of  two  atmospheric  paths.  The  slant 
path  originating  at  an  altitude  of  20  km  and  extending  to  space  at  a 75-deg 
zenith  angle  is  typical  of  viewing  geometries  for  satellite  surveillance 
systems."  The  100-km-long  horizontal  path  at  20-km  altitude  is  represen- 
tative of  the  viewing  geometry  used  in  on-going  field  measurements  pro- 
grams employing  sensor-equipped  aircraft.  An  equivalent  isothermal 
H2O/CO2  cylin’rical  plume  model  is  used  as  the  20-km  altitude  target. 
Results  for  this  source  over  the  two  atmospheric  paths  indicate  that  the 
effective  transmittance  of  the  atmosphere  to  hot  gas  emission  sources  may 
be  several  times  smaller  than  transmittance  values  that  are  calculated 
ignoring  the  line  correlation  effects. 
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I.  INTRODUCTION 


Consideration  of  the  infrared  transmittance  of  the  atmosphere  is  of 
fundamental  importance  to  the  analysis  of  the  spectral  characteristics  of 
missile  and  aircraft  plumes  viewed  through  long  atmospheric  paths.  For  ex- 
ample, in  order  to  design  surveillance  system  ser,so:3,  a £riori  knowl- 
edge is  required  of  the  expected  radiance  that  would  be  presented  to  such 
sensors.  Similarly,  methods  are  required  for  correcting  the  apparent 
radiance  measured  by  operating  systems  to  arrive  at  the  unattenuated  target 
signals.  The  solution  of  both  of  these  problems  hinges  on  the  dlity  to 
compute  appropriate  atmospheric  transmittance  factors. 

Quasi-empirical  methods  developed  for  computing  atmospheric 
attenuation  within  the  framework  of  band  model  approximations  have  been 
used  for  a number  of  years.  An  extensive  review  of  many  of  these  models 
is  given  by  Andmg.  Although  the  exact  procedure  used  in  developing  these 
various  models  differs  from  model  to  model,  the  underlying  approach  is  the 
same.  Basically,  a transmittance  function  T[k(v),  u]  is  determined  for  a 
particular  atmospheric  absorbing  species  from  laboratory  data  measured 
for  homogeneous  absorbing  paths  of  this  same  species.  The  variable  u is 
a measure  of  the  optical  depth  of  the  path,  and  k(v)  is  a measure  of  the 
mean  absorption  in  an  interval  Av  about  the  spectral  position  v.  Either  a 
theoretically  based  functional  dependence  of  t on  u and  k(v)  is  assumed  and 
k(v)  is  adjusted  for  a best  fit  with  the  laboratory  data  within  each  Aw  region, 
or  the  dependence  actually  displayed  by  the  data  is  determined  and  tabulated. 
With  either  approach,  th^  end  result  is  a transmittance  model,  which  can 
reasonably  predict  the  transmittance  spectrum  that  would  be  experimentally 


D.  Anding,  Band- Model  Methods  for  Computing  Atmospheric  Slant- Path 
Molecular  Absorption,  714Z-21-T,  Willow  Run  Laboratories.  Ann  Arbor. 
Mich.  (February  1967)(reprinted  January  1969). 


obtained  for  homogeneous  paths  of  arbitrary  optical  thickness.  In  order  to 
apply  these  models  to  general  atmospheric  slant  paths,  a means  is  required 
of  reducing  the  nonhomogeneous  slant  path  absorbing  properties  to  those  of 
an  equivalent  homogeneous  path.  This  reduction  is  usually  performed  by 
an  application  of  the  Curtis-Godson  approximation.  The  resulting  equivalent 
homogeneous  properties  are  then  inserted  into  the  band  model  to  yield  the 
transmittance  of  the  atmospheric  path. 

Anding*  has  made  a detailed  comparison  of  the  predictions  of  the  most 
prominent  of  these  band  models  and  the  limited  amount  of  actual  field  measure- 
ments that  exist  for  atmospheric  slant  paths.  In  general,  when  appropriate 
band  models  are  chosen  with  regard  to  specific  species  and  spectral  region, 
a reasonably  accurate  representation  of  atmospheric  transmittance  can  be 
obtained. 

Despite  the  reasonable  accuracy  of  these  models,  however,  two 
severe  limitations  preclude  their  application  to  the  present  problem.  First, 
the  models  are  applicable  only  for  optical  paths  over  which  the  temperature 
is  near  atmospheric  temperatures,  and  second,  they  are  appropriate  only 
for  predicting  the  attenuation  of  continuum  radiation  (e.g.,  blackbody 
radiation).  For  a great  number  of  targets  of  interest,  the  source  radiation 
is  not  continuum  radiation,  but  is,  in  fact,  highly  structured  (spectrally 
speaking)  molecular  band  emission.  The  temperature  limitation  imposed  on 
existing  models  prevents  their  application  to  the  computation  of  transmission/ 
absorption  properties  within  the  hot  gaseous  plumes  of  a missile  or  jet 
aircraft.  Moreover,  the  limitation  to  continuum  radiation  sources  prevents 
the  use  of  these  models  in  computing  the  transmission  of  radiation  emitted 
from  a particular  species  in  the  plume  when  this  same  species  is  present  as 
an  atmospheric  absorber  (e.g.  , H20).  The  high  degree  of  correlation 
between  the  line  positions  of  the  hot  gaseous  emission  spectrum  and  the  cool 
atmospheric  absorption  spectrum  results  in  a greater  amount  of  absorption 
than  would  be  effected  for  continuum  source  radiation. 
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The  inapplicability  of  existing  atmospheric  transmittance  models  for 

computing  hot  gaseous  radiation  characteristics  is  not  a serious  problem 

since  other  transmittance  models  explicitly  designed  for  hot  combustion 

2 

gas  calculations  have  been  formulated.  However,  the  inapplicability  of 
current  atmospheric  transmittance  models  to  calculations  involving  non- 
continuum source  radiation  is  not  so  easily  remedied. 

One  of  the  fundamental  objectives  of  the  present  work  is  to  assess 
the  degree  of  error  inherent  in  the  use  of  present  atmospheric  transmittance 
calculations  when  they  are  applied  to  the  problem  of  deattenuating  radiance 
data  obtained  for  targets  whose  primary  radiation  is  due  to  hot  gas  emission. 
At  present,  this  work  is  concerned  exclusively  with  missile  and  jet  aircraft 
plumes  containing  hot  CO^  and  H^O  as  the  only  infrared-active  species. 
Similarly,  the  only  atmospheric  attenuating  species  considered  are  CO^ 
and  cons*deration  given  to  particulate  emission  within  the  plume 

or  particulate  extinction  (absorption  or  scattering)  in  the  atmosphere.  The 
primary  combustion  products  for  hydrocarbon-fueled  missiles  and  aircraft 
are  H^O  and  CO.,,  and  these  are  the  radiating  species  for  which  the  line 
correlation  effect  would  render  the  most  inaccurate  calculations  using  present 
transmittance  models. 

The  line  correlation  effect  will  be  most  significant  for  the  attenuation 
of  hot  H^O  emission  by  cool  atmospheric  absorption  and  for  the  attenua- 

tion of  hot  CO^  emission  by  cool  atmospheric  CC^  absorption.  Since  the 
line  positions  between  H^O  and  CO^  are  uncorrelated,  however,  this  effect 
is  not  present  for  the  attenuation  of  hot  H^O  emission  by  atmospheric  CO^ 
absorption  or  vice  versa.  For  these  latter  conditions,  the  atmospheric 
transmittance  computed  using  present  band  models  would,  in  principle,  be 
correct.  Similarly,  present  band  models  would  be  applicable  for  calculating 
the  atmospheric  attenuation  of  radiation  emitted  from  particulate  matter  in 

! 

4 

2 

C.  B.  Ludwig,  W.  Malkmus,  J.  E.  Reardon,  and  J.  A.  L.  Thompson, 
Handbook  of  Infrared  Radiation  from  Combustion  Gases,  NASA  SP-3080, 
ed.  R.  Goulard  and  J.  A.  L.  Thompson  (1973). 
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plains s fe  g.  , carbon)  cr  from  “.lgine  hot  parts  since  this  radiation  would 
be  continuum  in  nature.  It  should  be  emphasized  that  this  line  correlation 
effect  exists  only  within  the  concept  of  a band  model  formulation,  i.e., 
for  the  case  where  the  average  transmittance  effects  in  a nonzero  spectral 
region  Av  containing  man/  individual  vibration-rotation  lines  arc  considered. 
Cn  a strictly  monochromatic  basis,  the  transmittance  of  the  atmosphere  is 
a function  of  the  atmospheric  properties  alone  and  is  not  influenced  by  the 
source  spectral  characteristics. 


The  present  work  is  concerned  with  the  design  and  application  of  a 
band  modeling  technique  that  includes  the  effects  of  line  correlation 
between  the  emitting  and  absorbing  species.  This  is  accomplished  by  treat- 
ing the  entire  line  of  sight  extending  from  the  senuor  position  through  the 
intervening  aunosphere  to  the  far  boundary  of  the  hot  gas  emission  source 
as  a single  radiating/absorfcing  path.  Thus,  the  atmospheric  and  source 
paths  are  treated  in  a coupled  manner.  The  effective  transmittance  of  the 
atmosphere  to  the  line  emission  source  is  derived  by  dividing  the  radiance 
computed  for  this  total  path  by  the  radiance  computed  for  the  source  portion 
of  the  path  alone  This  procedure  will  always  yield  a transmittance  factor 
smaller  than  the  atmospheric  transmittance  computed  for  a continuum 
source.  The  penalty  for  treating  the  problem  in  this  manner  is  the  genera- 
tion cd  a nigh ly-  irmomogeneous  line  of  sight.  Although  the  atmospheric  path 
pcrt'.or!  cf  the  total  path  is  itself  inhomogeneous,  the  degree  of  inhomcgenei  ly 
eai,  oe  .'Adequately  candled  witnin  the  Curtis- Godson  approximation.  Coupling 
x hot  gat  emission  source,  to  the  end  of  Ihe  atmospheric  path,  however,  can 
generate  a high  degioe  cf  inhomoger.ei'y.  The  principal  region  of  discon- 
tinuity occurs  at  the  boundary  of  t.ie  atmospheric  and  source  portions  of  the 
path.  As  will  be  shown  in  this  discussion,  the  degree  of  inhomogeneity 
generated  in  this  region  can  easily  exceed  the  limitations  of  the  Curtis- 
Godfjon  approximation.  To  circumvent  the  failings  of  the  Curtis-Godson 
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approximation  in  these  areas,  a recently  developed  formulation,  the 
Lindquist-Simmons  approximation,  is  invoked.  Results  derived  using  the 
latter  approximation  show  significant  improvement  over  those  obtained 
using  the  Curtis-Godson  approximation. 

Section  II  of  this  report  is  devoted  to  the  formulation  of  the  band 
model  used  in  the  present  work  and  to  a detailed  investigation  of  bofn  the 
Curtis-Godson  and  T indquist-Simmons  approximations.  A discussion  of 
the  band  model  parameters  used  within  the  model  is  presented  in  Section  III, 
and  a brief  account  of  the  computer  code  ATLES  (Atmospheric  Transmittance 
to  Line  Emission  Sources)  used  to  perform  the  transmittance  calculations 
is  given  in  Section  IV.  Results  obtained  for  hypothetical  missile  plume 
sources  for  two  typical  viewing  geometries  are  presented  in  Section  V.  A 
summary  and  recommendations  for  future  work  are  given  in  Section  VI. 


G.  H.  Lindquist  and  F.  S.  Simmons,  "A  Band  Model  Formulation  for  Very 
Nonuniform  Paths,"  J.  Quant.  Spectrosc.  Radiat.  Trans.  12,  807  (1972). 


II.  JV  ND  MODEL  FORMULATION 


A.  RADIANCE  EQUATIONS 

The  exact  solution  for  the  radiance  at  position  s = 0 of  a general  line 
of  sight  extending  to  c = S at  a specific  wave  number  and  for  a single 
optically  active  gas  is 


5 

L(v)  = / L*( v,  s)  k(v,  s)  c(s)  p(s) 

•A 


-f  k(v,  s')  c(s')  p(s')  ds' 


ds 


a) 


where  L*(v,  s)  is  the  Planck  radiation  function,  k(v,  s)  is  the  spectral 
absorption  coefficient  of  the  active  gas,  p(s)  is  the  total  gas  pressure,  c(s) 
is  the  mole  fraction  of  the  active  gas,  s is  the  geometric  path  length  variable, 
and  S is  tile  total  geomi  +ric  path  length.  This  formulation  of  the  equation 
assumes  that  there  is  no  source  of  radiation  (background)  beyond  s = S. 

The  variables  entering  into  this  equation  are  written  explicitly  as  functions 
of  the  path  variable  s to  emphasize  the  fact  that  they  may  change  with  path 
position. 

The  monochromatic  treatment  of  radiance  with  the  use  of  Eq.  (1)  for 
an  entire  vibration- rotation  band  of  an  active  molecule  might  require 
radiance  calculations  for  several  thousand  individual  lines.  In  addition,  a 
number  of  calculations  would  be  required  over  the  profile  of  each  line  to 
account  for  its  spectral  shape.  If  high  spectral  resolution  results  are  not 
required,  a considerable  reduction  in  calculation  effort  can  be  effected  by 
considering  the  average  radiation  effects  in  a small  interval  Av  containing  a 
band  of  several  tens  or  hundreds  of  lines. 

The  average  radiance  in  a spectral  interval  Av  centered  around  v is 


L(v) 


&L  L,v,,dv' 


(2) 
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Preceding  page  blank 


IfAv  is  small  enough  so  that  L*(w,  s)  is  sensibly  constant  over  Aw,  Eqs.  (1) 
and  (2)  can  be  made  to  yif  "d 


L(W>  = "f  L*<v>  s)  "Tjs>  ^ dt 


(3) 


where  t (w,  s)  is  the  average  transmittance  in  At 


(W.S)  = — f T(V',S)  d 

Av 


and  t(v,s)  is  the  spectral  transmittance 


t(w,  s)  = e 


/ k(w,  s/)c(8')  p(s/)  dj 


(4) 


(5) 


Equation  (3)  is  the  basis  of  radiance  calculations  using  band  model 
approximations.  Rather  than  attempting  to  characterize  the  radjation  trans- 
fer for  all  wave  numbers  of  interest  through  the  use  of  Eq.  (1),  band  model 
approximations  attempt  to  characterize  the  average  transfer  effects  in  a 
small  interval  Aw.  The  main  effort  in  band  modeling  is  the  characterization 
of  the  transmittance  term  t(w,  s)  entering  into  Eq.  (3)  without,  of  course, 
resorting  to  its  calculation  via  Eqs.  (4)  and  (5). 

The  generalization  of  Eq.  (3)  to  the  case  of  several  active  species  is 
made  by  writing  the  average  total  transmittance  7 as  a product  of  the  indi- 
vidual average  transmittances  of  each  species 


— m _ 

t(v,s)  = /7  t.  (w,s) 
i = 1 1 


(6) 
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where  m is  the  number  of  species.  This  relation  requires  the  assumption 
that  there  exists  a complete  lack  of  correlation  between  the  line  strengths 
and,  particularly,  the  line  positions  of  any  two  species  in  the  interval  Av. 
This  lack  of  correlation  has  been  demonstrated  to  be  true  for  H^O  and  CO^ 

Substitution  of  Eq.  (6)  into  Eq.  (3)  yields  the  following  form  of  the 
radiance  equation  used  in  the  present  work: 


/S  _ m 

L*(v,  s)  t(v,  s)  Y]  -=^7 — 
u iTl  Ti[v’ 


s) 


d s) 
ds 


ds 


(7) 


To  compute  the  transmittance  of  an  atmospheric  path  to  a line- 
emission  source,  radiance  calculations  using  Eq.  (7)  are  made  over  two 
paths.  The  first  path  is  the  source  path  alone,  and  the  second  path  is  the 
entire  line  of  sight  consisting  of  the  atmosphere  and  the  source  path  com- 
bined. Calculation  over  the  first  path  yields  the  radiance  of  the  source  at 
the  source  boundary,  while  calculation  over  the  second  path  yields  the  total 
radiance  at  the  sensor  (s  = 0)  position,  taking  into  account  the  absorption/ 
emission  effects  presented  by  the  intervening  atmosphere.  The  true  trans- 
mittance of  the  atmosphere  to  the  source  is  computed  by  dividing  the  total 
path  radiance  by  the  source  path  radiance.  This  true  transmittance  can 
then  be  compared  to  the  transmittance  t given  by  Eq.  (6).  The  latter 
transmittance  is  the  normally  computed  atmospheric  transmittance  and  is 
appropriate  to  a source  that  emits  continuum  radiation  (e.g.,  a blackbody) 
with  approximately  constant  intensity  across  Av.  The  effect  of  a high  degree 
of  correlation  between  the  line  positions  of  the  emitting  source  and  those  of 
the  absorbing  atmosphere  results  in  the  true  transmittance  being  less  than 
ti  e transmittance  to  a constant  continuum  source. 


4 

R.  M.  Goody,  Atmospheric  Radiation,  I:  Theoretical  Dasis,  Clarendon 
Press,  Oxford  (1964). 
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B.  BAND  MODEL 


Equation  (7)  is  the  relation  used  to  compute  radiance  in  the  band  model 
approximation.  The  aim  of  the  band  model  proper  is  to  characterize  the 
average  transmittance  t to  be  used  in  this  relation.  The  present  work  is 
carried  out  within  the  framework  of  the  t .atistical  band  model.  In  this 
model,  the  arrangement  of  lines  (i.e.,  the  line  positions)  within  a spectral 
interval  is  assumed  to  be  completely  random.  Within  this  approximation, 
the  average  transmittance  for  a single  species  is  given  by 


/ W\  f i N 

V T)=exP‘N6E  W(i) 

i=  1 


where  W is  the  mean  equivalent  width  for  all  the  lines  in  the  spectral  region 
Av,  6 is  the  mean  line  spacing  in  Aw,  and  N is  the  number  of  lines  in  Aw. 

To  provide  a working  model  when  the  individual  W(i)  are  not  all  known,  some 
analytic  intensity  distribution  is  generally  assumed  and  applied  to  Eq.  (8). 
Commonly  used  distribution  functions  include:  (1)  a constant  line  intensity 
distribution  wherein  all  lines  in  Aw  are  assumed  to  have  the  same  strength, 

(2)  an  exponential  line  intensity  distribution,  and  (3)  an  inverse  line- strength 
distribution  function.  For  reasons  to  be  discussed  later,  the  present  work 
has  been  carded  out  for  a random  array  of  pressure-broadened  Lorentzian 
lines  of  constant  line  intensity  When  this  distribution  is  assumed,  the  mean 

transmittance  for  a random  array  of  pressure-broadened  Lorentzian  lines  is 
given  by 

' * ,9, 

where  f(x)  is  the  Ladenburg-Reiche  function,  x is  the  dimensionless  optical 
depth  parameter,  and  (3  is  the  mean  line  width/ spacing  band  model  parameter. 
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A d'scussion  of  x and  p is  presented  in  paragraph  II.  C.  An  equivalent  way 
to  write  Eq.  (9)  is 


T'  - 1 

t = e 


(10) 


where 


T1  = 1 - pf(x) 


(11) 


The  term  t'  is  interpreted  to  represent  the  mean  transmittance  in  Av  for  a 
band  of  randomly  arranged  but  well- isolated  lines.  The  function,  .orm  of 
Eq.  (10)  then  allows  one  to  write  the  transmittance  for  general  overlapping 
conditions  in  terms  of  t'.  Note  that,  in  this  interpretation,  it  is  entirely 
possible  for  7'  to  be  a negative  number.  This  occurs  when  the  condition  of 
nonoverlapping  of  lines  no  longer  holds,  and,  under  these  conditions,  it  is 
best  not  to  think  of  7'  as  a transmittance,  but  simply  as  an  intermediate 
function  needed  to  compute  7.  Note  also  that  negative  values  of  t*  are,  in 
fact,  required  to  allow  t to  assume  values  of  less  than  e — 0.  368. 

In  addition  to  7,  the  derivative  d-r/ds  is  required  in  the  application  of 
Eq.  (7).  From  Eqs.  (10  and  (11), 


dT  _ — dr' 
ds  ds 


= -p  ^ - f (x)  ^ 


ds 


dx  ds 


ds 


(12) 


(13) 


C.  INHOMOGENEOUS  PATH  CONSIDERATIONS 

For  a homogeneous  path,  the  dimensionless  optical  depth  parameter 
x is  given  by 


x<8>  = y (i4) 

where  u = cps  is  the  optical  depth  of  the  path,  and  k and  0 ave  the  basic  band 
model  parameters.  The  term  k is  the  mean  absorption  coefficient  for  the 
interval  Av,  and  (3  is  a measure  of  the  ratio  of  the  mean  line  width  to  the 
mean  line  spacing  in  Av 


P=2"|  (15) 

For  an  inhomogeneous  path,  appropriate  path-averaged  values  for  k (or  x) 
and  p must  be  used.  In  the  present  work,  two  methods  of  accounting  for 
path  inhomogeneities  are  used:  In  addition  to  the  standard  Curtis -Godson 
approximation,  the  recently  developed  Lindquis t-Simmons  approximation 
is  considered.  The  latter  approximation  was  developed  to  correct  some  of 
the  failings  of  the  Curtis-Godson  approximation  when  applied  to  highly 
inhomogeneous  paths.  In  either  of  these  approximations,  the  following  path- 
averaged  quantities  are  used  for  x and  P: 

k(s')  c(s')  p(s')  ds'  (16) 


x (s)  = — 
e p 


7m  / 


0e<3>  = 


/ 


k(s')  c(s')  p(s')  p(s')  ds 


// 


k(s')  c(s')  p(s')  ds' 


(17) 


In  the  Curtis-Godson  approximation,  these  path-averaged  values  arc  used 
directly  for  x and  p in  Eqs.  (11)  and  (13)  to  give  the  nonoverlapping  trans- 
mittance and  transmittance  derivative 


t'  = 1 - p f(x  ) 
e e 


(18) 
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(19) 


dr'  . df(xe)  ** 

37  = ”Pe  “37 31 

e 


“37  ■ f(xe)  in 


Differentiation  of  Eqs.  (16)  and  (17)  allows  Eq.  (19)  to  be  written  in  the 
convenient  form 


= - k cp  y(xe,p) 


where 


y(xfP)  = (2-p)  + (p-  l) 


and  where  p = P/ p is  the  ratio  of  the  local  p to  the  path-averaged  value  up 
e 

to  that  point. 

In  effect,  the  Curtis- Godson  approximation  simply  replaces  the 
inhomogeneous  path  by  an  equivalent  homogeneous  path  whose  :c  and  P vari- 
ables are  the  path-averaged  variables  x^  and  Pe>  The  Lindqui:  t- Simmons 
approximation  also  uses  the  path-averaged  variables  x^  and  p^  but  in  a 
lesser  degree  of  approximation.  Rather  than  making  approximAtions  in  the 
expression  for  the  transmittance,  the  Lindquist- Simmon 5 approach  i3  to 
make  approximations  in  the  expression  for  the  transmittance  deri,rativc. 

The  rationale  for  this  approach  is  that  the  radiance  equations,  i.  e.,  Eq.  (3) 
or  (7),  are  formulated  basically  in  terms  of  the  transmittance  derivative 
rather  than  the  transmittance  itself.  Thus,  an  approximation  to  the  former 
should  be  inherently  more  accurate  than  an  approximation  to  the  latter.  For 
a single,  isolated  Lorentzian  line,  the  transmittance  derivative  is 


c(s)  p(s) 


-«  Y (s)  + Av 


1 rs  S(s')  Y(s') 
exp / — r 

”J0  Y^s')  + Aw 


c(s')  p(s')  ds'  d(Av>) 
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where  S and  y are  the  line  strength  and  width,  respectively.  For  a band  of 
nonoverlapping  Lorentzian  lines  of  equal  strength  and  width,  the  correspond- 
ing band  model  average  transmittance  can  be  written 


dT' 

nr  = 


i r k(8)  3(s) 

" -«  p2(s)  + Z2 


c(s)  p(s ) exp 


i re  k(s-)  p(S>) 

p2(s')  + Z2 


c(s')  p(s')  ds' 


dZ 


(23) 


where  Z = 2trAv/6.  The  necessity  of  assuming  a band  of  equal- strength 
lines  to  arrive  at  Eq.  (23)  is  one  of  the  reasons  that  the  constant  line  inten- 
sity distribution  form  of  the  statistical  band  model  has  been  assumed  through- 
out. The  need  for  this  restriction  will  be  the  subject  of  future  work.  In  the 
Curtis-Godson  approximation,  the  path-dependent  parameter  p(s)  occurring 
in  Eq.  (23)  is  replaced  by  the  path-averaged  value  Pfi.  This  has  the  effect  of 
performing  a separation  between  the  path  integration  and  spectral  integration 
operations  and  results,  after  performing  the  spectral  integration,  in 
Eqs.  (18)  and  (19).  In  the  Lindquist-Simmons  approximation,  the  replace- 
ment of  p by  pe  is  made  only  in  the  exponential  term  of  Eq.  (23),  which  can 
then  be  simplified  to 


dT' 

ds 


k cp  y(xe,p) 


(24) 


where 


y(x,  P)  = & f 


-x(  1 + cos  G) 


* J0  (p2  + 1)  + (p2  - 1)  cos  0 


dG 


(25) 


where,  G - 2 tan  * (2irAv/Pe6)  and,  again,  p = P/Pe>  Thus,  a complete 
separation  of  spatial  and  spectral  integrations  is  not  effected  in  the  Lindquist- 
Simmons  approximation.  A detailed  analysis  of  the  function  y is  given  in 
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paragraph  II.  D.  Note  also  that,  in  the  Lindquist-Simmons  approximation, 
no  analytical  solution  is  given  for  t;  and,  hence,  this  function  must  be 
integrated  numerically  after  dr'/de  is  evaluated  along  the  path,  i.e.. 


:'(e)  = 1 - f 

Jo 


<fr'  (s) 


A summary  of  formulas  required  for  implementing  the  band  model 
calculations  in  both  the  Curtis-Godson  and  Lindquist- Simmons  approximations 
is  /iven  in  Table  1.  The  Ladenburg-Reiche  function  f(x)  is  approximated  by 


f(x)  = 


1 + 


df(x) 

dx 


-4  (?r]  [<♦  (?)] 


-(1  + 1/ 2a) 


where  a = 5/4.  This  approximation  has  been  stated  to  be  accurate  to  less 
than  1 percent  for  all  x.  5 

The  Planck  radiation  function  used  in  Eq.  (7)  is 


L*( v,  s)  = 1.  192  X 10"  12v3  [e1,  439v/T(8)  _*] 


and  gives  L*(v,  s)  in  units  of  W/cm  sr  cm"  1 for  v in  units  of  cm"  1 and  T(s) 
in  units  of  °K. 


A.  Goldman,  "On  Simple  Approximations  to  the  Equivalent  Width  of  a 
Lorentz  Line,"  J.  Quant.  Spectrosc.  Radiat.  Trans.  8,  829  (1968). 
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Table  1.  Summary  of  Band  Model  Computation  Formulas 


Table  1.  Summary  of  Band  Model  Computation  Formulas  (Continued) 


Path  Averages 


y-Function 


Ladenburg-  Reiche 
Function 


Planck  Function 


D.  ANALYSIS  OF  y(x,  p) 

The  essential  difference  between  the  standard  Curtis-Godson  (CG) 
and  Lindquist- Simmons  (LS)  approximations  is  the  formulation  for  the 
function  y.  A greater  degree  of  approximation  characterizes  the  CG  than 
the  LS  approximation  and  allows  for  a complete  separation  of  the  spatial 
and  spectral  variables  of  integration.  The  LS  approximation  does  not  yield 
a complete  separation  of  spatial  and  spectral  variables  and,  thus,  requires 
a new  spectral  integration  for  each  new  spatial  coordinate  considered  along 
the  inhomogeneous  path.  An  analysis  vf  the  function  y has  been  made  for 
both  the  CG  and  LS  approximations  with  the  intent  of  showing  under  what 
conditions  the  CG  approximation  can  be  expected  to  fail,  how  the  LS  approxi- 
mation succeeds  under  these  same  conditions,  and  how  the  calculation  of 
y in  the  LS  approximation  is  effected  in  practice. 

The  function  y is  dependent  on  the  two  parameters  x and  p (the  sub* 
script  cn  xg  is  suppressed).  The  term  x is  a measure  of  the  optical  depth 
from  the  beginning  of  a line  of  sight  to  a point  along  the  line  of  sight,  and  p 
is  the  ratio  of  the  local  band  model  parameter  p to  the  average  value  of  p 
over  the  path  up  tc  that  point.  The  parameter  x is  zero  at  the  beginning 
(sensor  position)  of  any  optical  line  of  sight  and  generally  becomes  larger 
as  the  distance  along  the  line  of  sight  increases.  There  are  conditions, 
however,  under  which  x can  decrease  as  the  distance  into  the  optical  path 
increases.  If  the  defining  expression  for  x(s)  in  Eq.  (16)  is  differentiated 
with  respect  to  the  path  variable  s,  it  can  be  derived  that 

[2  - p(.H  ,30) 

Thus,  at  any  point  where  p is  greater  than  2,  dx(s)/ds  is  negative  and  x(s) 
decreases  with  increasing  distance  into  the  path.  The  largest  value  of  x 
encountered  in  any  calculation  to  date  has  been~2000. 
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The  parameter  p has  the  value  of  unity  at  the  beginning  of  any  line  of 
sight,  and  the  variation  of  p along  the  line  of  sight  then  depends  on  the  path 
type  and  the  strength  of  the  inhomogeneiiies  along  the  path.  For  atmospheric 
paths,  two  conditions  prevail.  If  a line  of  sight  originates  at  a high  altitude 
and  progresses  deeper  into  the  atmosphere,  p increases  gradually  and 
generally  approaches  a final  value  of  ~2.  However,  if  a line  of  sight 
originates  low  in  the  atmosphere  and  progresses  upward  to  low  pressure 
regions,  p may  obtain  values  as  small  as  10"  . For  a line  of  sight  that  pro- 
ceeds from  a cool  region  into  a hot  region,  p may  increase  substantially 
near  the  boundary  of  the  region  and  result  in  a local  increase  in  p up  to 
values  of  100  or  more.  The  largest  value  encountered  in  calculations  to 
date  has  been ‘'-200.  Thus,  it  is  evident  that  there  is  a large  range  over 
which  the  parameters  x and  p can  vary  and  this  range  has  important  impli- 
cations in  the  calculation  of  y,  particularly  within  the  LS  approximation. 

Regardless  of  the  extent  of  the  variation  of  x and  p,  however,  certain 
physical  restraints  can  be  placed  on  the  functional  form  of  a realistic  y- 
function.  These  restraints  are:  y — 1 as  x— 0 (weak  line  or  linear  limit); 
dy/dx— • Oasx— *0;y~x  as  x — co  (strong  line  or  square  root  limit); 
y — ► 1 as  p -*a>;  and  y remains  less  than  or  equal  to  unity  for  all  x and  p. 

1.  CURTIS- GODSON  APPROXIMATION 

In  the  CG  approximation,  y(x,  p)  is  given  by  Eq.  (21)  where  f(x)  is 
approximated  by  Eq.  (27).  From  these  two  equations, 


Differentiation  of  this  function  with  respect  to  x and  solving  for  the  extrema 
yield  the  conclusion  that  y will  have  an  extremum  (maximum)  for  x > 0 
whenever  p 2 2 + 1 fnr  = 2.80.  The  position  of  the  maximum  is 
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'max 


«(p  - 2) 


i + 


| [•!<>  - 2)  - >]} 


■i  + 1/2  a 


(33) 


Since  y(x  = 0)  is  unity,  y must  therefore  become  greater  than  unity  for 
x > 0 whenever  p > 2 + i /a.  The  position  at  which  y again  becomes  less  than 
unity  can  be  obtained  by  solving  the  equation  y(x,  p)  = 1.  One  solution  is 
obviously  x = 0,  and  the  other  solution  is  contained  in  the  expression 


El*\  = 


" + *?)-( 


1 + 


irx\ 

T) 


3/2 


(34) 


For  p > 10,  this  condition  requires  x > p /2ir  before  y i3  again  less  than 
unity.  Some  curves  for  y(x,p)  in  the  Curtis-Godson  approximation  are 
shown  in  Fig.  1 . 

The  failing  of  the  CG  approximation  can  thus  be  seen  to  occur  primarily 
under  conditions  where  a strong  inhomogeneity  gives  rise  to  a large  increase 
in  p (large  p)  near  tire  beginning  of  a line  of  sight  (small  or  medium  x).  This 
ic  the  condition  that  obtains,  for  example,  for  a calculation  on  a line  oi 
sight  originating  on  the  boundary  of  a rarified  missile  plume  or  flame.  If  no 
large  increase  in  p occurs  until  far  into  the  path  where  x is  large,  no  inaccu- 
racy is  immediately  apparent  in  the  CG  approximation.  However,  under 
these  conditions,  the  region  of  large  inhomogeneity  contributes  little  to  the 
line -of- sight  radiance  because  of  the  small  transmittance  of  the  previous 
portions  of  the  path.  The  result  of  Eq.  (30)  is  important  here  also.  If  p is 
large  enough  to  make  the  CG  approximation  inaccurate,  it  is  also  large 


-24- 


enough  to  make  dx/ds  negative.  Thus,  even  with  increasing  depth  along  the 
geometric  path,  x does  not  increase  and  the  path  calculations  become  trapped 
in  the  region  of  maximum  error. 

2.  LINDQUIST-SIMMONS  APPROXIMATION 

A distinct  advantage  of  the  CG  approximation  is  that  y is  expressible 
in  terms  of  f(x)  and  the  latter  function  can  be  approximated  by  a rather 
simple,  but  quite  accurate,  analytical  expression.  The  y-function  in  the 
LS  approximation,  however,  cannot  be  written  so  simply.  The  calculation 
of  y in  the  LS  approximation  by  a new  numerical  spectral  integration  for 
each  path  interval  position  is  highly  inefficient  if  accuracy  of  better  than 
approximately  1 percent  is  required.  The  approach  to  the  use  of  the  LS 
appro:rimation  hae  been  to  compute  y numerically  over  a grid  of  x,  p values, 
and  then  to  perform  a two-dimensional,  log-log,  second-order  interpolation 
on  the  table  to  obtain  y at  specific  x,  p points.  In  addition,  analytic  approxi- 
mations have  been  found  for  certain  conditions  of  x and  p . 

The  exact  expression  for  y(x,  p)  in  the  LS  approximation  is 


(x,p  ) = ^/  7— 

J0  (1  + P ) + 


-x(l  + cos  0) 


f0  (1  + p ) + (1  - p ) COS  e 


No  analytic  solution  for  this  integral  has  been  discovered.  An  analysis  of 
the  integrand  of  this  integral  shows  that,  in  addition  to  the  extrema  that 
always  occur  at  0 = 0 and  0 = it,  the  integrand  has  an  extremum  between 
0 = 0 and  0 = tt  for  p < 1 only  and  then  only  for 


i 


-26- 


This  extremum  is  a minimum  and  only  one  extremum  exists.  These 
conditions  delineate  the  three  cases  for  the  positive  x,  p quadrant  shown  in 
Fig.  2(a).  The  form  of  the  integrand  between  6 = 0 and  6 = it  for  each  case 
is  shown  in  Fig.  2(b).  Note  that  the  condition  of  very  large  x in  cases  I 
and  III  leads  to  the  formation  of  a very  sharply  peaked  function  around  0 = tt, 
while  the  condition  of  very  small  p in  cases  II  and  III  leads  to  a very  sharply 
peaked  function  around  0 = 0.  These  conditions  result  in  the  ability  to  find 
analytic  asymptotic  solutions  for  the  extreme  conditions  of  these  cases. 

The  asymptotic  solution  for  very  large  x has  been  obtained  using  a standard 
analysis  procedure,  which  results  in  the  correct  square  root  form 


y(x, p) — 


J 


2ir 


+ 1p  jlR 
+ 2 


I 


(37) 


The  conditions  for  accuracy  of  this  asymptotic  expression  have  been  deter- 
mined by  comparison  to  numerical  calculations  and  are  given  in  Table  2. 


Table  2.  Conditions  for  Accuracy  of  the  Large-x  Asymptotic  Solution 
for  y(x,  p)  in  the  Lindquist-Simmons  Approximation 


% Accuracy 

Minimum  Allowed  x 

0.01 

x £ 1300 

. 2 
and  x 4 85  p 

0.  10 

130 

8.7  p2 

1.0 

14 

2. 5 p2 

10. 

5.  3 

0. 33  p2 

For  very  large  p also. 

the  integrand  is 

sharply  peaked  about  0 = it  . 

However,  the  asymptotic  solution  obtained  for  large  p is  neither  realistic 
nor  accurate,  and  the  solution  is,  in  fact,  the  same  as  that  for  large  x. 
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(b)  FORMS  OF  THE  INTEGRAND 

Fig.  2.  Integrand  of  y(x,  p)  Integral  in  the 
Lindquist -Simmons  Approximation 
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For  large  p and  small  x,  ideally  y should  approach  unity  whereas  the 
solution  obtained  approaches  1 This  effect  is  probably  due  to  the  fact 

that  the  integrand  is  not  a strongly  enough  varying  function  with  p around 
0 = -rr  to  render  the  asymptotic  technique  valid. 

A similar  condition  was  found  when  attempting  to  extract  an  asymptotic 
solution  for  small  p.  In  this  case,  a sharply  peaked  function  occurs  around 
0=0.  The  asymptotic  solution  obtained  is 


y 


(38) 


Although  this  result  is  independent  of  p,  it  is  unrealistic  in  that  it  also 

approaches  1 fyfv  rather  than  unity  as  x— »-0.  However,  it  was  discovered 

/—  . 2x 

that,  by  neglecting  the  factor  1/yir,  the  remaining  function,  namely,  e , 

provides  an  extremely  good  approximation  to  y for  small  p and  small  x. 

This  approximation  can  be  obtained  also  by  approximating  the  numerator  of 

the  integrand  for  small  x and  small  p by  e , i.e.  , its  value  at  0 = 0 where 

most  of  the  contribution  to  the  integral  is  concentrated.  For  small  p and 

large  x,  a sharp  peak  occurs  around  both  6=9  and  9 = n (case  III),  and  an 

accurate  solution  can  be  obtained  by  adding  the  two  asymptotic  solutions  for 

these  two  conditions  together. 


y — » e"2x  + P.  — » asp-0  (39) 

FMS 

The  upper  limit  on  p,  for  which  this  solution  is  accurate  to  a given  percent- 
age error,  ha3  not  been  considered  in  detail.  A comparison  with  ex»  •. 

_4 

numerical  calculations  at  p = 10  , however,  h?s  shown  the  maximum  error 

to  be  of  the  order  of  1 percent  with  this  maximum  error  occurring  at  x ^ 5. 
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For  p = 1C"4,  this  value  of  x is  the  region  where  the  two  contributions  to  the 
asymptotic  solution  ere  of  the  same  magnitude  (Fig.  3)  Attempts  to  adjust 
the  asymptotic  solution  obtained  for  large  p by  factors  of  tt  or  /n  applied 
at  various  places  in  the  equation  did  not  lead  to  any  improved  form  for  an 
asymptotic  solution. 

For  10"4  < x < 104  and  10"4  < p < 104,  numerical  calculations  of 
y(x,p)  were  performed  and  the  results  used  to  form  an  interpolation  table 
as  discussed  earlier.  Calculations  were  performed  for  both  x and  p at  the 
25  values:  1 X 10"4,  2 X 10"4,  5 X 10"4,  1 X 10  3,  2 X 10  3,  ^ 

5 X 103,  and  104.  The  accuracy  of  the  calculations  was  set  at  1 part  in  10 
although  some  entries  may  be  more  accurate.  The  upper  limit  on  x and  the 
upper  and  lower  limits  on  p were  set  by  the  fact  that  the  number  of  A9 
intervals  required  for  accurate  numerical  integration  increased  extremely 
xapidly  beyond  these  limits.  The  numerical  integration  was  performed  using 
the  Wedlc  approximation  wivh  the  number  of  intervals  between  0 and  ir  equal 
to  1 2n  where  n is  an  integer.  The  limiting  value  of  n < 200  set  the  x and  p 
limits  given  above.  The  calculated  y values  are  listed  in  Table  3 and  illus- 
trated in  Fife.  3.  A summary  of  the  calculatlonsl  procedure  for  all  x and  p 
is  given  in  Fig.  4. 

A comparison  of  the  y-functions  appropriate  to  the  CG  and  LS 
approximations  reveals  the  following* 

a.  In  both  cases,  y(x,p)  is  well  behaved  for  very  large  x and 
app^o-schfes  the  correct  asymptotic  form  y — p/./  2t\x. 

h.  The  functions  v(x,  p)  are  identical  for  the  case  p = 1 (as 

they  should  be  since  this  value  reflects  the  conditions  for 
^ homogeneous  path). 

c.  In  the  CG  approximation,  y can  exceed  unity  for  path 
conditions  where  p xl  2 + l/o.  In  the  LS  approximation, 
y never  exceeds  unity. 

d.  The  transition  from  the  small- x behavior  to  the  asymptotic 
large-x  behavior  generally  occurs  very  much  faster  in  the 
U3  approximation  than  in  the  CG  approximation.  Except  for 
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very  large  p,  x has  achieved  the  asymptotic  form  between 
* ~ 3 ”*.8  *n  LS  approximation.  However,  for  small  p 
in  the  CG  approximation,  this  changeover  to  the  asymptotic 
form  may  be  prolonged  as  far  as  x—  102  or  103  or  even 
further. 
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abulation  of  y(x,  p)  in  Lindquist-Simmons  Approximation  for  x and  p 


Tabic  3.  Tabulation  of  /(x,  p)  in  Lindquiat-Simmona  Approxim**’,.,!  for  > : >"i 
Bstwron  0. 0001  and  10,000  (Continued) 


Table  3.  Tabulation  of  y(x,  p)  in  Lindquist-Simmons  Approximation  for  x and  p 
Between  0.0001  and  10,000  (Concluded) 
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Fig.  4.  Summary  of  y(x,  p)  Calculation  in  the 
Lindquist -Simmons  Approximation 


in.  BAND  MODEL  PARAMETERS 


The  basic  data  required  for  performing  band  model  calculations  are 
the  band  model  parameters  k and  (3  (or,  equivalently  for  (3,  the  average  line 
width  parameter  V and  the  mean  line  spacing  parameter  6 since  P = 2ttY/6). 
These  parameters  must,  of  course,  be  consistent  with  the  band  model 
formulation  used.  Three  sets  of  these  parameters  are  presently  in  use. 

A.  GENERAL  DYNAMICS  PARAMETERS 

A set  of  band  model  parameters  that  was  generated  by  General 

Dynamics  and  intended  primarily  for  use  in  high-temperature  applications 

has  been  in  existence  for  several  years.  The  most  recent  publication  of 

2 

these  parameters  is  given  by  Ludwig  et  al.  For  both  CO^  and  H^O,  the 
listed  variables  are  k and  6 and  art  appropriate  for  use  in  the  statistical 
band  model  formulation  with  an  exponential  line  intensity  distribution. 

The  data  for  CO^  are  the  result  of  theoretical  calculations  based  on 
the  known  quantum  mechanical  properties  of  linear  triatomic  molecules. 
Results  obtained  using  these  data  are  generally  conceded  to  be  in  good  agree- 
ment with  experimental  results.  The  k parameter  for  CO^  is  complete  for 
both  the  4.3-fj.m  and  2.7-p.m  bands  for  temperatures  between  300  and  3000°K. 
The  6 data  for  the  4.  3-p.m  band  are  also  complete  for  the  temperature  range 
of  300  to  3000°K.  However,  for  the  2.  7-fim  band,  6 data  are  given  only  up 
to  1800°K.  Both  the  k and  6 data  are  given  at  spectral  intervals  of  5 or 
10  cm"*  and  presumably  reflect  values  appropriate  to  spectral  resohition 
intervals  of  about  the  same  size. 

General  Dynamics  data  for  H^O  are  derived  almost  exclusively  from 
experimental  measurements.  Because  of  the  nonlinearity  of  H^O,  the  exact 
theoretical  calculation  of  band  model  parameters  is  extremely  complicated 
and  difficult,  parti cularly  at  high  temperatures.  The  bulk  of  the  data 
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above  ~1200°K  were  obtained  in  a consistent  manner  from  emission/ 

absorption  measurements  made  on  very  long  strip-burner  H^/O^,  flames. 

Data  below  ~1200°K  are  based  on  extrapolations  from  these  high-temperature 

data  and  on  the  analysis  of  published  low-temperature  H,0  spectra. 

- 1 z 

Recommended  values  are  given  in  25  cm  steps  over  the  spectral  region 
from  50  to  9300  cm"  * for  k and  1150  to  7500  cm  * for  6.  The  spectral 
resolution  width  is  also  25  cm’*.  The  k data  are  given  for  the  temperature 
range  from  300  to  3000°K,  but  the  6 data  are  given  from  only  600  to  3000°K. 

The  generation  of  p data  from  the  General  Dynamics  6 data  is  effected 
using  the  line  width  parameters  given  in  Table  4 (which  has  been  taken  from 
Ref.  2).  The  line  width  (half- width  at  half-height)  for  the  species  i is 
computed  from  these  parameters  by 


Y.(T) 


i'. »,  m"1  j 


■ <Vo  p. 


(40) 


The  summation  term  accounts  for  foreign  gas  broadening  (j  * i)  and  non- 
resonant self-broadening  (j  = i)  contributions  to  the  line  width.  The  second 
term  accounts  for  the  resonant  self-broadening  effects  which  occur  for  H^O. 
The  term  p.  is  the  partial  pressure  of  the  ith  species.  Data  for  k,  6,  and  Y 
have  been  used  to  generate  a card  file  of  k(v,  T)  and  P_(v,  T)  data  for  H,G 

1 U L 

and  CO,.  The  units  of  k are  cm  /atm  and  were  obtained  bv  multiplying  the 
c _ 1 
General  Dynamics  data  by  273/T.  The  units  of  P^  are  atm  for  nonresonant 

self-broadening  only.  Simple  semilogarithmic  linear  extrapolations  were 

used  to  generate  6 values  for  H^O  at  temperatures  below  600°K  and  Tor  CC^ 

at  temperatures  both  above  1800°K  and  below  300°K. 

Although  the  General  Dynamics  data  are  presumably  appropriate  for 
the  statistical  band  model  with  the  assumption  of  an  exponential  intensity  dis- 
tribution, these  data  are  actually  quite  suitable  for  use  with  the  assumption 
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Table  4.  General  Dynamics  Line  Width  Parameters 


Molecule  (i) 

Broadener  (j) 

(Vi>0 

(cm"  * / atm) 

n. 

J 

* 

>0 

(cm"  * /atm) 

h2o 

H2° 

0.  09 

1/2 

0.44 

C°2 

0.  12 

1/2 

N2 

0.09 

1/2 

°2 

0.04 

i/2 

C°2 

C°2 

0.09 

1/2 

0 

H2° 

0.07 

1/2 

N2 

0.07 

1/2 

°2 

0.  06 

1/2 

VT)  * £ (<*{>„  Pj  {W)  J}  + <\>oPi(- 

J 

* 

n. 

273\  1 
T) 

V 

j 

> 

Nonresonant  self- 
broadening and 
foreign  gas 
broadening 


Resonant  self- 
broadening 


of  a constant  line  intensity  distribution.  The  conventional  curve  of  growth 
for  the  exponential  distribution  function  is 


W 

6 


■ P 


x 


(41) 


In  order  for  this  function  to  assume  the  correct  form  in  the  limit  of  strong 
absorption,  the  parameter  (3  has  to  be  defined  as  4/it  times  the  parameter  (3 
consistent  for  use  with  the  constant  line  intensity  curve  of  growth.  The 
constant  line  intencity  p-parameter  can  be  used  with  the  exponential  distri- 
bution curve  of  growth  if  this  factor  4/n  is  absorbed  into  th-  functional 
dependence  of  W/6.  When  this  functional  transformation  is  performed,  the 
resulting  exponential  distribution  curve  of  growth  becomes 


(42) 


Since  this  is  the  curve  of  growth  function  recommended  for  use  with  the 
General  Dynamics  band  model  parameters,  it  is  concluded  that  the  General 
Dynamics  parameters  are  suitable  for  use  with  the  constant  line  intensity 
distribution  curve  of  growth  also. 

B.  LINE- AVERAGED  PARAMETERS 

The  second  and  third  cels  of  band  model  parameters  were  generated 
from  the  comprehensive  line  compilation  atlas  compiled  by  AFCRL.  ^ 

This  data  set  provides  line  data  for  each  of  over  100,000  lines  for  seven 
atmospheric  absorbing  species  inching  HzO  and  C02-  The  data  given  for 


R.  A.  McClatchey,  W.  S.  Benedict,  S.  A.  Clough,  D.  E.  Burch, 

B f * Calfee,  K.  Fox,  L.  S.  Rothman,  and  J.  S.  Garing,  AFCRL  Atmos- 
pheric Absorption  Line  Parameters  Compilation.  AFCRL-TR-7  3-006, 

Air  Force  Cambridge  Research  Laboratories,  Mass.  (25  January  1973). 
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each  line  include  the  line  position  v(cm”*)»  the  line  strength  S (cm  */ 
molecule  citi  at  296°K),  the  line  half- width  parameter  V (cm”  /atm  for 
air  broadening  at  296°K),  and  the  energy  of  the  lower  level  of  the  transition 
E (cm-1).  Two  approaches  have  been  taken  in  averaging  these  line  data 
over  small  spectral  intervals  to  provide  band  model  parameters.  The  first 
approach  is  the  straightforward  interpretation  that  the  appropriate  averages 
are  given  simply  by 


1 ■ srZ  s,i) 

i=  1 


00'  — 


where 


l **-’/ 


In  these  expressions,  Av  is  the  width  of  the  spectral  interval,  N is  the 
number  of  lines  in  the  interval,  and  6 is  the  mean  line  spacing  6 = Av/N. 

The  zero  subscript  on  p and  V indicate  values  per  unit  pressure.  The  second 
approach  provides  for  an  exact  fit  in  the  weak  line  and  strong  line  limits 
between  the  general  statistical  band  model  formulation  and  the  statistical 
band  model  in  which  the  constant  line  intensity  distribution  function  is 
assumed.  This  development  is  outlined  by  Goody.  In  the  general  statistical 
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band  model,  the  average  integrated  absorptance  for  a band  of  isolated 
spectral  lines  is 


<4t) 

i=l 

where  W(i)  is  the  equivalent  width  for  the  ith  line  in  the  inverval.  For  a 
Lorentzian  line  and  homogeneous  path, 

W(i)  = 2n  Y(i)  f(x)  (47> 

where  f(x)  is  the  Ladenburg-Reiche  function,  x is  the  dimensionless  optical 
depth  parameter  x = S(i)  u/2ir\(i),  and  u is  the  optical  depth.  Thus,  in  the 
general  statistical  band  model  formulation,  A is  given  explicitly  in  terms  of 
the  line  data  S(i),  Y(i).  When  the  constant  line  intensity  distribution  is 
applied  to  the  statistical  band  model,  the  average  integrated  absorptance 
appropriate  to  the  interval  Av  for  a band  of  isolated  spectral  lines  is 

X = P f(x)  <48) 


where 


(49) 


In  this  case,  the  average  integrated  absorptance  is  given  in  terms  of  the 
band  model  parameters  k and  p.  In  the  weak  line  limit,  as  x — 0 and 
f(x)  — x,  the  equating  of  A given  by  Eqs.  (46)  and  (48)  then  results  in  an 
expression  for  k that  is  the  same  as  Eq.  (43).  In  the  strong  line  limit, 
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as  x— »oo  and  f(x)  — * y/Zx/ir,  the  equating  of  A given  by  Eqs.  (46)  and  (48) 
leads  to  the  result  given  by  Eq.  (44)  except  that  '/  is  now  given  by  the 


expression 


y = JL  LllL 

T0  N 


E./ SilTyi) 


N 

E S(i) 

i=  i 


The  essential  difference  between  these  two  approaches  is  that  the 
second  approach  applies  a line- strength  dependent  weighting  operation  to  the 
line  widths.  This  second  approach  is  significant  in  that  it  imparts  a tem- 
perature dependence  (through  the  temperature  dependence  of  S)  to  Yn  other 

-1/2  U 
than  just  the  normal  T dependence  that  is  applied  in  the  first  approach. 

These  line -averaging  calculations  are  performed  by  the  computer 
program  LINAVE,  which  reads  in  the  AFCRL  line  data  and  performs  the 
averaging  calculations  in  specified  spectral  intervals  over  a specified  tem- 
perature range.  The  line  strengths  S(i)  are  computed  from  the  value 
given  at  the  standard  temperature  Tq  = 296  K by 


S.(T)  = S.(TQ) 


Qr  <to>  qv(to) 
qr  (T)  Qv(T) 


-v(i)/kT 
- v(i)/kTr 


lil/l  . J_ 
\T  To 


The  factors  QD  and  Q.r  are  the  rotational  and  vibrational  partition  functions 
respectively.  The  exponential  term  accounts  for  the  Boltzmann  distribution 
population  of  molecules  in  the  energy  level  E(i),  and  the  term  in  brackets 
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accounts  for  stimulated  emission  effects.  The  ratio  of  rotational  partition 
functions  is  given  by 


qr  <to> 

o^nr 


(52) 


where  m = 1 for  H^O  (H^O  is  assumed  to  be  a nearly  symmetric-top 
molecule)  and  m = 3/2  for  CO^.  The  molecular  vibrational  partition  func- 
tion is  assumed  to  be  the  product  of  simple  harmonic  oscillator  partition 
functions  for  each  vibrational  mode  of  the  molecule. 


“di  -d  -d_ 

r -v  /kTl  r -v  /kTl  - v /kTl  3 

QV(T)  = [l-e  J [l-e  J j^l  - e J (53) 

whore  v^,  and  are  the  fundamental  oscillatory  frequencies  of  the 
molecule  and  d^,  d^»  and  d^  are  the  degeneracies  associated  with  these 
modes.  For  CO^t  = 1388.  17  cm  \ vz  = 66’7.40  cm"\  and 

v3  = 2349.  t6  crn"1.  For  H20,  = 3657.0  cm"1,  vz  = 1594.7  cm*1,  and 

v3  = 3755.7  cm  . The  degeneracies  for  all  modes  of  both  CO  and  H O 
are  unity,  except  for  the  doubly  degenerate  bending  mode  in  the  linear 
co2  fur  which  d2  = 2.  The  value  1/k  = 1. 439  cm  °K  is  used 

throughout. 

Tl-.o  temperature  variation  given  in  Eq.  (51)  accounte  for  the  change 

i:i  the  population  of  absorbers  in  the  energy  level  E(i)  with  temperature. 

Another  temperature  variation  must  be  included  to  convert  S.  from  the 

• 1 — 2 i 

AFCRL  units  of  cm  /(molecule  cm  ) to  the  standard  units  used  in  this 
- 2 

work,  i.e.,  cm  /atm.  This  conversion  is  effected  by  multiplying  S.(T) 
obtained  from  Eq.  (51)  by  the  conversion  factor 


2.  480  X 1019 


molecules 
atm  cm3 


(54) 
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Tho  line  broadening  parameters  included  in  the  AFCRL  line  atlas 
are,  presumably,  the  broadening  parameters  appropriate  to  foreign  gas 
(air)  broadening.  The  line  width  data  of  Table  4 have  been  used  to  convert 
these  parameters  to  reflect  values  appropriate  to  nonresonant  self- 
broadening. Note  that  only  the  ratios  of  the  line  width  parameters  from 


Table  4 are  used  and  not  the  line  width  parameters  themselves.  The  VQ(i) 


values  used  in  Eqs.  (45)  and  (50)  arc  thus  determined  from  the  AFCRL 
values  by  multiplying  the  latter  by  the  factor 


Z71 

c \ I~~y~ 


(55) 


where  c = 1.  186  for  H O and  c = 1.  382  for  CO  . 

^ 2 


The  LINAVE  Program  computes  k (v,T)  and  PQ  (v,  T)  where  k has 
units  (cm  /atm)  and  pQ  has  units  (atm-1  for  nonresonant  self-broadening). 
Most  of  the  work  presently  being  done  uses  k and  PQ  values  computed  for 
Av  = 25  cm  that  have  been  computed  for  wave  number  steps  of  5 cm’1 
between  2500  cm’  and  4400  cm’  1 for  HzO  and  between  3075  cm"  1 and 
3875  cm  for  CC>2  for  14  temperatures  between  100°K  and  3000°K. 

C.  COMPARISON  OF  PARAMETER  SETS 

The  k value  for  C02  at  the  2.  7-^m  band  center  position  v = 3700  cm’  1 
is  shown  in  Fig.  5.  In  this  comparison,  the  General  Dynamics  set  and  the 
line-averaged  set  are  seen  to  agree  well  up  to  a temperature  of  ~1000°K. 

A comparison  of  PQ  parameters  for  the  same  condition  is  shown  in  Fig.  6 
where  good  agreement  between  the  General  Dynamics  parameter  and  the 
line-averaged  parameter  computed  from  Eq.  (50)  is  evident  up  to  a tempera- 
ture of -600  K.  However,  a substantial  difference  may  be  seen  to  occur 
between  the  line-averaged  pQ  values  computed  from  Eqs.  (45)  and  (50).  The 
discrepancy  between  k and  pQ  from  Eq.  (50)  for  temperatures  higher  than 
-600  K can  probably  be  explained  by  the  fact  that  the  AFCRL  line  compilation 
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R icrrf'/atmi 


I 


LINAVE1 


was  generated  primarily  for  calculations  involving  atmospheric  slant  paths 
where  maximum  temperatures  of  only  ~350°K  are  likely  to  be  encountered. 
Therefore,  lines  that  become  important  only  at  higher  temperatures  have 
not  been  included  in  the  compilation. 

Similar  comparisons  for  H_0  at  wave  numbers  v = 3600  cm*'  (near 

^ - i 

the  2.  7-fjm  band  center)  and  v = 3200  cm  (in  the  2.  7 -pm  band  wing)  are 

shown  in  Figs.  7 and  8.  The  values  of  k from  the  two  sets  for  the  band 

center  region  (Fig.  7)  are  also  seen  to  compare  well.  In  the  band  wing 

position,  however,  there  is  little  resemblance  between  the  temperature 

variations  of  k from  the  two  sets.  For  temperatures  less  than  ~1000°K, 

the  most  probable  cause  of  the  discrepancy  in  the  wing  is  the  experimental 

error  involved  in  attempting  to  measure  very  small  absorptance  values. 

For  the  higher  temperatures,  the  most  probable  cause  of  the  discrepancy 

is  again  the  lack  of  inclusion  of  high-temperature  lines  in  the  AFCRL  line 

compilation.  These  same  arguments  can  also  be  applied  to  the  P_  compari- 

- i 0 - 1 

son  shown  in  Fig.  8(a)  for  H^O  in  the  3200  cm  region.  In  the  3600  cm 
regior\,  illustrated  in  Fig.  8(b),  a reasonably  good  transition  from  the  low- 
ternperature  line-averaged  p^  computed  from  Eq.  (50)  to  the  high-trmperature 
Gtneral  Dynamics  PQ  is  evident  in  the  region  of  1250°K.  A detailed  investi- 
gation of  the  matching  of  the  low-temperature  line-ave-aged  parameters  and 
the  high-temperature  empirical  General  Dynamics  parameters  for  H^C  will 
be  the  subject  of  future  work.  At  present,  one  or  the  other  of  these 
parameter  sets  is  used  at  all  temperatures. 

D PRO/ DENING  EFFECTS 

The  basic  band  model  parameter  data  input  used  in  the  fransmittance 
calculations  are  the  k and  p^  data  just  discussed.  The  p^  parameters  are 
the  intrinsic  values  for  nonresonant  self-broadening  with  units  of  atm"*. 

The  actual  value  of  the  band  model  parameter  P depends  on  the  particular 
species  composition  existing  at  the  point  where  P is  calculated.  For  both 
atmospheric  paths  and  hot  gas  source  paths,  the  only  species  considered  to 
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be  present  are  CO^*  H^O,  and  air  (78  percent  and  22  percent  by 
pressure).  From  Eq.  (40)  and  the  ratio  of  line  width  parameters  given  in 
Table  4,  the  parameter  p can  be  computed  in  terms  of  the  nonresonant 
self-broadening  parameter  p^.  For  H^O  and  CO^,  the  resulting  equations 
are,  respectively: 


P>2o(T) 


H-O 

= Pq  (T)  p 


[• 


0.  878  + 0.456  c 


CO. 


+ 0.  122  + 4.  89 


Mh2o 


(56) 


Pco<T> 


CO 

= pn  ^ (T)  p 10.753  + 0.  247  c 


0. 


CO. 


+ 0.0244  c 


»2°] 


(57) 


where  p is  the  total  pressure  in  atmospheres  and  an^  are 

mole  fractions  of  CO^  and  h^O  respectively.  For  atmospheric  paths,  only 
the  first  term  in  the  brackets  is  significant  in  each  case.  These  terms  are 
simply  the  inverse  of  the  c-factors  appearing  in  Eq.  (55). 
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IV.  COMPUTER  CODE  (ATLES) 


A computer  program  has  been  generated  to  handle  the  computation  of 
atmospheric  transmittance  to  line  emission  sources.  This  code  is  called 
ATLES  (Atmospheric  Transmittance  to  Line  Emission  Sources)  and  is 
written  in  FORTRAN  IV  compatible  with  The  Aerospace  Corporation 
CDC  7600  computer  and  operating  system.  A simplified  flow  diagram  for 
this  program  is  shown  in  Fig.  9.  The  following  listing  is  a brief  discussion 
of  each  of  the  numbered  blocks  on  the  flow  chart. 

1.  Each  run  of  the  program  requires  a set  of  data  cards  to 

control  the  calculational  mode  of  the  program.  Pertinent 

data  entered  in  this  set  include: 

a.  The  number  of  equal  length  intervals  that  the 
atmospheric  path  (NAPATH)  and  source  path 
(NSPATH)  are  divided  into  for  numerical 
integration  over  these  paths. 

b.  The  type  of  source  being  viewed.  (Present 
capability  provides  for  a constant  continuum  source, 
a single  line -of- sight  path,  and  an  axisymmetric 
cylindrical  plume  source.  For  the  latter  source, 
the  number  of  intervals  is  also  entered  into  which 
the  axial  and  transverse  coordinates  of  the  projected 
area  of  the  plume  are  divided  for  numerical 
integration  over  the  plume.) 

c.  The  viewing  geometry  of  the  sensor / source  configu- 
ration. (These  data  include  sensor  and  source 
altitudes,  zenith  angles,  and  ranges  as  required.) 

A spherical  earth  geometry  is  used,  but  atmospheric 
refraction  effects  are  not  included. 

d.  The  spectral  data  for  the  minimum  (vrnin)  and 
maximum  (v'max'  wave  number  and  the  number  of 
spectral  intervals  (NSPEC)  to  be  considered. 

e.  A series  of  parameters  that  direct  the  program  to 
read  in  new  large  blocks  of  data  or  to  repeat  the 
data  used  in  tha  previous  run  (if  there  was  one),  anu 
a set  of  parameters  that  direct  the  program  to  print 
out  or  suppress  the  printout  at  intermediate 
calculation  stages. 
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Fig.  9.  Flow  Diagram  of  ATLES  Frogram  (Continued) 
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Fig.  9.  Flow  Diagram  of  ATLES  Program  (Concluded) 
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2.  The  data  defining  the  model  atmosphere  to  be  used  in  the 
calculations  are  re'id  in  by  subroutine  MODATM.  The 
format  for  these  data  provides  for  the  entering  of 
atmospheric  pressure,  temperature,  H2O  concentration, 
and  aerosol  concentration  at  altitudes  of  0,  1,  2,  3,  . . . 

24,  25,  30,  35,  . . . 50,  70,  and  100  km.  This  altitude 
format  is  consistent  with  the  model  atmospheres 
presented  by  McClatchey  et  al.,7  that  have  been  used 
for  all  computations.  The  McClatchey  data  are  converted 
to  units  consistent  with  the  operational  units  of  the  pro- 
gram. Pressure  is  converted  from  millibars  to 
atmospheres  via 

platm)  = P(miuib*r> 

1.013  X 10J 


H2O  vapor  concentration  is  converted  from  density  in 
gm/m^  to  mole  fraction  via 

ph  o T (°K> 

Cjj  q = 4.559  * 10”^  p(atm) 


CO2  concentration  is  set  at  cqq2  = 0.00033. 

The  data  entered  at  altitudes  specified  above  are  used  to 
construct  an  altitude  interpolation  table  with  the  thermo- 
dynamic variables  given  at  each  1-km  step  between  0 
and  100  km. 

3.  Subroutine  VIEW  uses  the  viewing  geometry  entered  in 
Step  1 to  compute  the  altitude  variation  along  the  line  of 
sight  from  the  sensor  to  the  source  at  NAPATH  + 1 
equally  spaced  points  along  the  path.  Portions  of  a path 
above  100-km  altitude  (or  an  entered  limit  ZMAX)  are 
not  considered. 


' 




R.  A.  McClatchey,  R.  W.  Fenn,  J.  E.  A.  Selby,  F.  E.  Volz,  and 
J.  S.  Caring,  Optical  Properties  of  the  Atmosphere  (Revised), 
AFCRL-71-0279,  Air  Force  Cambridge  Research  Laboratories,  Mass. 
(10  May  1971). 
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Subroutine  APATH  uses  the  altitudes  derived  in  Step  3 
and  the  model  atmosphere  interpolation  table  set  up  in 
Step  1 to  compute  pressure,  temperature,  and  concen- 
trations (pTc  data)  at  the  NAPATH  + 1 points  along  the 
atmospheric  path. 

This  step  is  skipped  if  calculations  for  a constant 
continuum  source  are  being  made.  If  a single  line  of 
sight  (LOS)  is  being  considered,  the  pTc  data  at 
selected  positions  along  the  source  path  are  read  in  by 
subroutine  PLUME  1.  If  a cylindrical  plume  is  being 
considered,  tho  pTc  data  at  selected  radial/axial  grid 
points  are  read  in  by  subroutine  PLUME2.  These 
data  constitute  interpolation  arrays  for  subsequent 
calculations . 

The  senscr  spectral  response  function  is  read  in  by 
subroutine  SENSOR. 

The  bard  model  parameters  k and  Pq  for  both  H^O  and 
CO2  a^e  read  in  by  subroutine  BAND.  This  set  con- 
sists of  a wave  number  array,  a_temperature  array, 
and  a two-dimensional  array  of  k and  pQ  for  the  tem- 
perature and  wave  numbers  of  the  first  two  arrays. 
These  data  also  constitute  interpolation  tables  for 
subsequent  calculations. 

A header  card  on  this  data  block  identifies  the 
parameters  so  that  conversion  routines  can  be  applied 
to  the  data  to  provide  a set  of  parameters  with  con- 
sistent units.  The  standard  unit  for  k is  cm”  V atm  and 
the  standard  unit  for  Pq  is  atm*'  for  nonresonant 
s elf-broadening . 

The  vmin,  vmax  and  NSPEC  data  entered  in  Step  1 are 
used  to  initiate  a spectral  loop  over  the  NSPEC  + 1 
equally  spaced  wave  numbers  between  ^rnin  and  ^max- 
If  NSPEC  = 0,  calculations  at  vrnin  only  are  performed. 

From  the  temperature  variation  over  the  atmospheric 
path  computed  in  Step  4 and  the  wave  number  set  by 
Step  8,  subroutine  KBPATH  computes  k(^o?»  Pq^Oz, 
ku  q,  and  ovcr  the  atmospheric  path  by  inter- 

polation on  the  data  tables  generated  in  Step  7.  These 
Pq  values  are  further  processed  using  the  pTc  data 
from  Step  4 in  subroutine  BETA  to  yield  the  final  band 
model  parameter  p at  each  point  along  the  atmospheric 
path. 
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10. 


Using  the  k(3  values  generated  in  Step  9 and  with  the  pTc 
data  of  Step  4,  subroutine  DELTRAN  computes  the 
transmittance  tj  and  the  transmittance  derivative  dT^/ds 
along  the  atmospheric  path  for  i = and  CO2  indi- 

vidually in  both  the  Curtis -God  son  and  Lindquist-Simmons 
approximations.  The  function  y(x,p)  required  in  the 
Lindquist-Simmons  approximation  is  computed  in  the 
function  subprogram  YLS.  The  transmittance  variables 
are  then  used  in  subroutine  RADNCE  to  compute  the 
radiance  of  the  atmospheric  path.  Subroutine  RADNCE 
also  returns  the  path  ending  values  of  the  t ->tal  trans- 
mittance TH2O  X “CO?  » which  is  the  atmospheric  trans- 
mittance appropriate  tor  a constant  continuum  source. 

11.  If  a constant  continuum  source  is  being  considered,  a 
skip  to  Step  19  now  takes  place.  Otherwise,  subroutine 
SGRID  is  called.  This  subroutine  calculates  the  appro- 
priate looping  parameters  required  to  effect  a two- 
dimensional  integration  loop  over  the  projected  area  of 
a cylindrical  source  plume.  For  each  LOS  passing 
through  the  plume,  Steps  12  through  17  are  performed. 

If  only  a single  LOS  is  being  considered,  the  loop  is  gone 
through  just  once. 

12.  Using  the  looping  parameters  from  Step  11  and  the  viewing 
geometry  data  read  in  at  Step  1 (principally,  the  viewing 
aspect  angle),  subroutine  SPATH2:  (a)  computes  the 
geometry  of  intersection  between  the  LOS  and  the  plume, 

(b)  computes  the  axial  and  radial  coordinates  (in  the 
coordinate  system  of  the  plume)  for  each  of  NSPATH  + 1 
equally  spaced  points  along  the  portion  of  the  LOS  within 
the  plume,  and  (c)  uses  these  axial/radial  parameters  to 
effect  an  interpolation  on  the  pTc  data  table  read  in  by 
subroutine  PLUME 2 (Step  5)  to  arrive  at  the  pTc  variation 
over  the  source  path. 

13.  This  step  computes  the  k(3  parameters  over  the  established 
source  path  in  the  same  manner  as  these  parameters  were 
evaluated  for  the  atmospheric  path  except  that  the  pTc  data 
from  Step  12  are  used. 

14.  Again,  this  step  is  similar  to  the  calculations  of  Step  10 
except  that  the  data  appropriate  to  the  source  pa*n  are 
used. 

15.  Step  14  computes  the  radiation  characteristics  over  the 
source  path  portion  of  the  LOS  only.  In  Step  15,  the 
radiation  properties  of  the  total  path  (i.  e.  , the  combined 
atmospheric  path  plus  source  path)  are  computed.  This 
step  is  essentially  the  same  as  Step  14  except  that  different 
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initial  conditions  at  the  source  path  starting  point  are 
imposed.  These  initial  conditions  are,  in  fact,  the  ending 
path  conditions  computed  for  the  atmospheric  path  in 
Step  10.  The  imposition  of  these  conditions  accounts  for 
the  absorption  effects  of  the  intervening  atmosphere. 

The  atmospheric  contribution  to  the  total  radiance  could  be 
accounted  for  by  adding  the  source  path  radiance  computed 
with  these  initial  conditions  to  the  atmospheric  path 
radiance  computed  in  Step  10.  In  practice,  the  atmospheric 
radiance  from  Step  10  is  not  included  in  the  total  path 
radiance.  The  atmospheric  radiance  calculation  accounts 
for  only  the  atmospheric  path  between  the  sensor  and  the 
source.  Thus,  inclusion  of  this  term  would  represent  an 
incomplete  account  of  the  atmospheric  radiance  since  the 
contributions  from  the  atmospheric  path  beyond  the  source 
(and  also  any  other  background  beyond  the  source)  would 
not  be  included.  The  complete  neglect  of  the  atmospheric 
radiance  contributions  to  the  total  path  radiance  makes  the 
calculated  results  analogous  to  (and  comparable  with) 
experimental  measurements  for  which  background  signalc 
have  been  subtracted. 

16.  The  true  atmospheric  transmittance  for  the  source  path 
defined  for  the  current  LOS  is  computed  by  dividing  the  total 
path  radiance  (minus  the  atmospheric  and  background 
contributions)  computed  in  Step  15  by  the  source  path 
radiance  computed  in  Step  14. 

17.  The  loop  over  all  LOS  through  the  projected  area  cf  the 
plume  is  continued.  Within  this  loop,  an  integration  over 
the  entire  projected  area  of  the  plume  is  effected  to  arrive 
at  the  total  source  and  sensor  radiance. 

18.  The  effective  atmospheric  transmittance  of  the  entire  source 
is  computed  by  dividing  the  total  source  radiance  by  the 
total  sensor  radiance. 

19.  The  loop  ovei  all  wave  numbers  between  v . and  v 

. “ , min  max 

is  continued. 


i 


I 


20. 


The  generated  spectral  transmittance  and  radiance  data 
are  averaged  over  the  sensor  response  function  read  in  by 
subroutine  SENSOR.  This  averaging  is  performed  by  sub- 
routine BNDPASS.  Bandpass  transmittances  for  line 
emission  sources  are  computed  by  dividing  bandpass 
averaged  radiances  and  not  by  averaging  the  computed 
spectral  transmittance.  The  bandpass  averaged  trans- 
mittance for  a continuum  source  i£  computed  by  averaging 
the  spectral  transmittance. 


i 

i 
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21.  A final  summary  of  the  various  spectral  radiance  and 

transmittance  calculations,  bandpass  averages,  and  input 
data  identification  codes  is  printed  out.  An  option  allows 
for  a punched  deck  output  for  use  in  plotting  routines. 

All  of  the  input  data  and  the  results  of  the  computations 
at  each  major  computation  step  can  be  printed  out  at 
option.  Such  printouts  are  indispensable  for  program 
checkout  and  detailed  analysis  of  selected  runs.  In  the 
normal  operating  mode,  however,  these  printouts  are 
suppressed  as  many  of  them  are  located  within  imbedded 
loops  and  would  involve  unmanageable  volumes  of  paper. 

As  presently  written,  this  code  should  be  considered  primarily  as  a 
research  tool.  It  was  written  to  enable  the  analysis  of  the  many  facets  of 
the  general  problem  of  transmittance  calculation.  The  general  calculational 
approach  is  geared  to  this  purpose  rather  than  to  the  efficient  calculation 
that  would  be  required  for  a production  analysis  tool,  although  this  code  is 
certainly  much  more  efficient  than  codes  performing  high-resolution  mono- 
chromatic calculations.  Because  of  this  approach,  the  code  is  generally  in 
a state  of  continual  modification,  but  only  minor  future  modifications  (e.g.  , 
input/output  formats)  are  envisioned.  Because  of  this,  a detailed  descrip- 
tion and  listing  of  the  code  have  not  been  included. 


V.  RESULTS  AND  DISCUSSIONS 


A.  HOMOGENEOUS  PATHS 

Comparisons  of  band  model  spectra  and  experimental  spectra  for 

homogeneous  hot  gas  optical  paths  were  made  in  order  to  assess  the  degree 

of  accuracy  with  which  the  LINAVE  parameters  (computed  from  the  AFCRL 

atmospheric  line  atlas)  could  predict  hot  gas  spectra.  Comparisons  were 

B 

made  with  an  experimental  CO^  spectrum  of  Burch  and  Gryvnak,  a CO^ 
spectrum  of  Simmons  et  al. , ^ and  an  H£0  spectrum  of  Simmons  et  al. 

The  experimental  conditions  appropriate  to  these  three  spectra  are  listed 
in  Table  5. 

Table  5.  Homogeneous  Hot  Gas  Spectra  Experimental  Conditions 


Reference 

Gas 

P 

(atm) 

E9 

L 

(cm) 

g 

Burch  and  Gryvnak 

pure  CO^ 

0.997 

1200 

7.75 

9 

Simmons  et  al. 

pure  CO^ 

2.03 

1146 

60 

Simmons  et  al. 

pure  H^O 

0.937 

1040 

60 

I 

8D.  E.  Burch  and  D.  A.  Gryvnak,  Infrared  Radiation  Emitted  by  Hot  Gases 
and  Its  Transmission  Through  Synthetic  Atmospheres,  U-1929,  Aeronutronic 
Division,  Ford  Motor  Company,  Newport  Beach,  Calif.  (31  October  1962). 

9F.  S.  Simmons,  H Y.  Yamada,  and  C.  B.  Arnold,  Measurement  of 
Temneratui  e Profiles  in  Hot  Gases  by  Emission-Absorption  Spectroscopy, 

NASA  CR- 72491 , National  Aeronautics  and  Space  Administration 
(April  1969). 

10F.  S.  Simmons,  C.  B.  Arnold,  and  D.  H.  Smith,  Studies  of  Infrared 
Radiative  Transfer  in  Hot  Gases,  I:  Spectral  Absorptance  Measurements 
in  the"2.7p  Bands,  4613-91-T,  Willow  Run  Laboratories  (August  1965). 
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Results  for  the  CO^  spectrum  of  Burch  and  Gryvnak  are  illustrated 
in  Fig.  10.  The  band  model  results  shown  in  Fig.  10(a)  were  computed  with 
LINAVE  parameters  generated  for  a resolution  of  25  cm‘*  and  at  wave 
number  steps  of  25  cm”  * . The  results  given  in  Fig.  10(b)  were  generated 
with  parameters  of  the  same  resolution  but  at  steps  of  only  5 cm  * . Com- 
parison of  the  experimental  result  with  the  band  model  result  computed  with 
the  General  Dynamics  band  model  parameters  has  been  discussed  by 
Ludwig  et  al.  Agreement  of  the  experimental  spectrum  with  those  com- 
puted using  the  LINAVE  parameters  is  surprisingly  good;  between  3550  and 
3800  cm"*  (the  band  center),  the  degree  of  agreement  is  about  as  good  as 
that  obtained  using  the  General  Dynamics  parameters.  In  the  wing  region 
below  approximately  3450  cm  *,  however,  the  results  for  the  LINAVE 

parameters  fall  seriously  short  of  the  experimental  curve.  Comparison 

q 

with  the  CO2  spectrum  of  Simmons  et  al.  (Fig.  1!)  displays  the  same 
features  just  discussed,  particularly  the  vast  underestimating  of  the  radiance 
in  the  wing  region  by  the  band  model  results  using  the  LINAVE  parameters. 
(Note  that  the  actual  experimental  curve  is  not  shown  in  Fig.  11  since  iv"  is 
negligibly  different  from  the  band  model  results  obtained  using  the  General 
Dynamics  parameters.)  The  most  probable  explanation,  as  mentioned 
earlier  in  paragraph  III.  C,  is  the  lack  of  inclusion  of  hot  line  data  in  the 
AFCRL  line  compilation. 

Comparisons  of  the  band  model  results  compwted  w;th  the  LINAVE2 
parameters  and  the  results  obtained  from  a monochromatic  line-by-line 
calculation**  for  these  two  CO  2 3pectra  are  presented  in  Figs.  10(b)  and 
1 1(b).  The  line-by-line  results  were  computed  on  a spectral  scale  of 
0.02  cm”*  and  then  convoluted  with  a 20-cm  * -wide  rectangular  slit  function 
for  comparison  with  the  band  model  results.  Line  data  required  for  the 
monochromatic  calculations  were  obtained  from  the  same  data  set  from 


C.  M.  Randall,  Monochromatic  Transmittance/ Radiance  Computations, 
TR- 007 5(5647)  - 3,  The  Aerospace  Corporation,  El  Segundo,  Calif. 

(31  December  1974). 
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(b)  Comparison  of  B«»d  Model  and  Line  by  Line  Results 


Fig.  iO.  Comparisons  of  Band  Model,  Line -by-Line,  and 
Experimental  Results  of  Burch  and  Gryvnak 
for  CC>2 
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which  the  LINAVE  band  model  parameters  were  derived.  Although 
agreement  between  the  band  model  and  line-by-line  spectra  is  quite  good, 
the  band  model  results  do  yield  a consistently  higher  degree  of  absorptance. 
Two  causes  may  be  responsible  for  this  effect.  First,  the  constant  line 
intensity  distribution  used  within  the  statistical  band  model  is  the  distribu- 
tion that  results  in  the  most  absorption.  If  the  intensity  distribution  were 

altered  to  allow  more  weak  lines  (e.g.  , an  exponential  distribution),  the 

2 

absorptance  would  become  less  and  agreement  with  the  line-by-line  results 

would  become  better.  Although  the  alteration  of  the  intensity  distribution 

to  allow  weak  lines  increases  the  number  of  strong  lines  also  (so  that  the 

mean  line  strength  remains  constant),  the  reduced  absorption  by  the  weak 

lines  outweighs  the  increased  absorption  of  the  strong  lines.  In  reference 

to  the  choice  of  band  model,  note  that  the  statistical  band  model,  regardless 

of  intensity  distribution,  is  generally  not  applied  to  CO_  absorption.  Most 

c 1 

often,  a regular  line  spacing  model  (e.g.,  the  Elsasser  model)  is  applied. 
However,  since  the  regular  model  does  not  allow  for  the  possibility  of  large 
gaps  between  lines  (as  is  allowed  in  the  statistical  model),  the  absorptance 
for  the  regular  model  is  greater  than  that  given  by  the  statistical  model. 

Thus,  the  use  of  a regular  line  spacing  model  would  lessen  the  agreement 
with  the  line-by-line  results. 

Secondly,  the  line-by-line  results  are  computed  in  a procedure  where- 
by the  exact  monochromatic  transmittance  is  first  computed  and  then  averaged 
(convolved  with  the  slit  function)  over  a spectral  region  Aw.  The  band  model 
approach,  on  the  other  hand,  first  performs  the  averaging  procedure  within 
Aw  (the  generation  of  the  band  model  parameters)  and  then  generates  the 
transmittance  spectrum.  There  is  no  reason  to  believe  that  the  two  mathe- 
matical operations  of  generating  the  spectrum  and  averaging  are 
commutative. 

The  comparison  of  the  band  model  calculations  using  the  LINAVE 

parameters  and  the  H-O  spectrum  of  Simmons  et  al.  is  shown  in  Fig.  12. 

C - 1 

The  band  model  parameters  were  calculated  for  a 25-cm  resolution  and 
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at  stops  of  25  cm’*.  Agreement  of  the  results  for  the  LINAVE2  parameters 
and  the  experimental  result  is  quite  good  above  about  3400  cm  How- 
ever, agreement  of  the  results  for  the  LINAVE1  parameters  is  less  marked. 
The  biggest  apparent  discrepancy  between  the  emissivity  results  computed 
using  the  LINAVE1  and  LINAVE2  parameters  occurs  for  those  conditions 
where  the  optical  path  is  not  thin  (c  ~ 0. 2)  but  is  not  large  enough  to  force 
the  emissivity  in  either  calculation  to  unity.  Since  most  of  this  H^O  s c- 
trum  is  not  thin,  the  LINAVE2  parameters  would  be  expected  to  yield  the 
better  results  since  these  parameters  were  derived  to  give  the  correct 
answer  in  the  strong  line  limit.  For  both  CO^  and  H^O,  the  LINAVE2 
parameters  are  used  throughout  the  remainder  of  the  calculations  presented 
in  this  report.  Results  of  the  band  model  calculation  using  a LINAVE2 
parameter  set  computed  for  a resolution  of  25  cm  , but  at  steps  of  only 
5 cm”*,  is  shown  in  Fig.  13  along  with  the  experimental  spectrum  and  the 
band  model  result  obtained  using  the  General  Dynamics  parameters. 

The  most  surprising  result  obtained  for  the  H^O  cornParison  is  the 
very  large  discrepancy  occurring  in  the  band  model  calculations  using  the 
LINAVE  parameters  for  the  wing  region  below  3300  cm  . This  result  can 
be  traceable  only  to  the  fundamental  line  data  from  which  the  LINAVE 
parameters  were  derived,  a conclusion  that  is  confirmed  by  the  nearly 
monochromatic  line-by-line  result  shown  in  Fig.  14.  This  spectrum  was 
computed  from  the  same  basic  line  data  from  which  the  LINAVE  parameters 
were  derived  and  also  shows  the  enhanced  features  in  the  wing  The  error 
in  the  line  data  tabulation  in  this  case  is  obviously  not  an  omission  of  lines 
(as  for  CO^)  but  either  an  inclusion  of  lines  chat  do  not  belong  or  an  error  in 
the  data  for  lines  that  do  occur  there.  Fortunately,  for  most  of  the  work  of 
concern  in  this  report,  only  the  spectral  region  above  3300  cm  * is  impor- 
tant. The  LINAVE2  parameters  can  thus  be  considered  to  be  reasonably 
accurate  up  to  at  least  1050°K. 
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Figure  14  contains  £.  comparison  of  the  band  model  calculations  and 
the  nearly  monochromatic  line-by-line  calculations  for  the  H^O  spectrum. 

As  in  the  case  for  CO^,  the  band  model  calculation  yields  a consistently 
higher  degree  of  absorptance  than  does  the  line-by-line  calculation.  The 
discussion  of  this  effect  given  for  CO.  applies  here  also.  A confirmation 
of  the  better  agreement  that  would  exist  for  a band  model  employing  a line 
intensity  distribution  that  allows  weak  lines  is  given  by  the  work  of  Goldman 
and  Kyle.  Agreement  to  better  than  5 percent  between  monochromatic 
line-by-line  calculations  and  band  model  calculations  employing  an 
exponential- tailed,  inverse  line  strength  distribution  was  achieved  by  these 
workers. 

B . ATMOSPHERIC  SLANT  PATHS 

Detailed  calculations  for  the  radiative  transfer  along  two  atmospheric 

paths  representative  of  geometries  appropriate  to  satellite- borne  and 

aircraft-borne  sensor  measurements  programs  were  performed.  The  path 

appropriate  for  a space  sensor  geometry  extends  from  the  source  altitude 

of  20  km  to  space  at  a zenith  angle  of  75  deg.  The  path  appropriate  to 

aircrafi-bornc  sensor  field  measurements  programs  consists  of  a 100-km- 

lorg  horizontal  path  situated  at  an  altitude  of  20  km.  The  hoi  gas  target 

source  used  fcr  these  calculations  was  based  on  a right  circular  cylinder, 

♦ 

horn  og  ere  ov.s  and  Isothermal  plume  model  given  by  Thomas.  The  plume 
dimensions,  temperature,  and  species  cc ncentraticns  for  plume  altitudes 
between  i and  63  km  are  listed  in  Table  6.  For  the  source  altitude  cf 
20  km,  L = 238  m,  r = 3.87  m,  T = 1250°K,  = 0.253,  and 

c Cq2  = 0.0371.  The  troprical  atmospheric  model  given  by  McClatchey  was 
used  in  these  calculations.  This  atmosphere  is  listed  in  Table  7.  Results 
were  calculated  using  both  the  LINAVE2  band  model  parameter  set 
(Av  = 25  cm  step  size  = 5 cm  S and  the  General  Dynamics  parameter  set. 

D.  D.  Thomas,  The  Aerospace  Corporation,  private  communication,  1974. 

^A.  Goldman  and  T.  G.  Kyle,  "A  Comparison  between  Statistical  Model  and 
Line-by-Line  Calculations  with  Application  to  the  9.  6y  Ozone  and  the  2.  7p 
Water  Vapor  Bands,"  Applied  Optics  7,  1167  (1968). 
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Table  6.  Homogeneous  and  Isothermal  Plume  Model 


z 

(km) 

r 

(m) 

L 

(m) 

P 

(atm) 

T 

(°K) 

1 

1.534 

20.  50 

0. 8924 

2139.  0 

5 

1.718 

37.64 

0. 5518 

1867.4 

10 

2.  170 

71. 35 

0. 2823 

1625.2 

15 

2.870 

133.9 

0.  1303 

1417.  8 

20 

3.865 

237.  9 

0. 05577 

1249.  4 

25 

5.197 

409.  5 

0.02537 

1107.  1 

30 

6.  897 

704.  1 

0.01204 

975.  8 

35 

8.918 

1. 260(3) 

5. 92 3 (—  3 ) 

844.4 

40 

11.29 

2.  237(3) 

3. 011(-3) 

723.5 

45 

14.09 

3.789(3) 

1 . 570(- 3 ) 

622.  5 

50 

17.54 

5. 944(3) 

8 430(-4) 

545. 6 

55 

22.46 

9.817(3) 

4.  574(-4) 

478.  7 

60 

29.41 

1.444(4) 

2.  382(-4) 

441. 5 

z = source  altitude 
r = plume  radius 
L = pi  ime  length 

p = pressure  (ambient  atmospheric  at  z) 
T = plume  temperature 


-j_j  q = mole  fraction 


CpQ  = C02  mole  fraction 
2 


-h2° 


CO, 
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Table  7.  Model  Tropical  Atmosphere 


Altitude 
(lan ) 

Pr  es  sure 
(atm) 

Temperature 

(°K) 

0 

1. 0000 

300.0 

1 

0 8924 

294.0 

2 

0.  *7947 

288.0 

3 

0.7053 

234.0 

4 

0.6249 

? 77.0 

i> 

0. 5518 

270.0 

6 

0. 4857 

264.  0 

7 

0.  4265 

257.0 

8 

0. 3731 

250.  0 

9 

0.  3248 

244.  0 

10 

0 2823 

237.0 

11 

0 Z438 

230.0 

12 

0 2103 

224.0 

13 

0.  179? 

217.0 

14 

0. 1540 

210.0 

15 

0. 1303 

204.0 

16 

0.  1096 

197.0 

17 

9.  250(-2) 

195.0 

18 

7 . 789(-?) 

199.0 

' o 

j 6.  575  (-  2) 

203. 0 

20 

5 577 (- 2) 

207.0 

21 

4.  738(*  2) 

211  0 

22 

4.  Q38(- i) 
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The  radiance  spectra  computed  for  the  homogeneous  source  alone  and 
for  the  sensor  position  are  plotted  in  Figs.  15  and  16  for  the  slant  path 
and  horizontal  path  respectively.  These  radiance  spectra  as  well  as  all  the 
results  presented  in  Figs.  15  through  24  were  computed  for  a single  line 
of  sight  that  intersects  the  equivalent  cylindrical  plume  at  a 90-deg  aspect 
angle  (measured  from  nose-on)  and  at  the  full  diameter  position.  Thus, 
the  source  is  a simple  homogeneous  path  2r  774  cm  in  length. 

For  the  LINAVE2  parameters,  the  sensor  radiance  spectra  computed 
using  the  CG  and  the  LS  approximations  are  nearly  identical.  The  reason 
for  this  result  is  that  the  LINAVE2  parameters  do  not  produce  a serious 
inhomogeneity  when  passing  from  the  cool  atmosphere  (~200°K)  into  the 
hot  source  (—1200  K).  It  may  be  seen  from  Fig.  6 that  for  CO?  at  3700  cm-*, 
p would  in  fact  obtain  a value  of  only  about  3 during  the  transition,  and  for 
H^O  (Fig.  8),  p changes  very  little  during  the  transition.  For  the  General 
Dynamics  parameters,  on  the  other  hand,  the  difference  in  the  results  when 
using  the  CG  approximation  or  the  LS  approximation  is  dramatically  evident. 
Here,  the  variation  of  the  band  model  parameter  p (see  Figs.  6 and  8)  with 
temperature  is  rapid  enough  to  provide  a very  strong  inhomogeneity  during 
the  transition  from  the  atmosphere  into  the  source  (p  ~ 100  for  CO^  and 
p - 5 for  H^O).  The  serious  breakdown  of  the  CG  approximation  is  evident 
in  its  prediction  of  a sensor  radiance  that  is  larger  than  the  source  radiance 
itself  at  some  spectral  positions.  This  physically  impossible  result  does 
not  occur  in  the  LS  approximation. 

Results  for  the  atmospheric  transmittance  t appropriate  for  a con- 
tinuum emission  source  are  presented  in  Figs.  17  and  18  for  the  slant  and 
horizontal  atmospheric  paths  respectively.  Comparisons  are  shown  of  the 
band  model  results  generated  using  the  LINAVE2  parameters  and  the 
General  Dynamics  parameters  as  well  as  results  computed  on  a monochro- 
matic line-by-line  basis**  and  results  obtained  using  the  LOWTRAN  band 
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rison  of  t Results  for  Horizontal  Path 
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model  routine.  In  general,  the  results  obtained  with  the  LINAVE2 
parameters,  the  line-by-line  calculations,  and  the  LCWTRAN  Program 
agree  quite  well.  The  causes  of  some  of  the  differences  that  do  occur 
between  the  LINAVE2  band  model  results  and  the  line-by-line  results,  par- 
ticularly in  the  regions  of  strongest  absorption,  are  probably  the  same  as 
those  responsible  for  the  discrepancies  in  the  homogeneous  path  spectri. 
discussed  in  paragraph  V.  A.  Another  possible  cause  lies  in  th ' fact 
that  the  band  model  results  were  performed  under  conditions  for  which 
necessary  numerical  integrations  along  the  optical  path  were  confirmed  to 
be  accurate  (see  paragraph  V.C).  However,  for  the  sake  of  reducing  com- 
puter processing  time,  the  line-by-line  calculations  were  performed  by 
dividing  the  atmosphere  into  only  5 to  10  equivalent  homogeneous  layers. 

This  number  of  layers  may  not  be  sufficient  to  obtain  high  accuracy  for  the 
path  integration  portions  of  the  calculations. 

Comparisons  of  t and  the  transmittance  appropriate  to  the  attenuation 

of  hot  gas  line  emission  t ^ are  shown  in  Figs.  10  through  22.  The  large 

difference  between  r and  caused  by  the  effect  of  line  correlation  between 

the  emission  and  absorption  spectra  is  evident  in  all  cases.  A transmittance 

ratio  of  nearly  4 is  displayed  at  3850  cm’  * for  the  horizontal  path.  The 

largest  differences  between  r and  generally  occur  in  the  regions  of 

strongest  absorption.  As  wa3  the  case  for  t,  the  difference  in  t between 

e 

calculations  involving  the  CC  and  15  approximations  is  very  small  when  the 
LINAVE2  parameters  are  employed.  When  the  General  Dynamics  parameters 
are  used,  however,  the  LS  approximation  yields  a decidedly  more  realistic 
result.  The  breakdown  of  the  Cvj  approximation  is  manifested  by  yielding 
effective  atmospheric  transmittance  values  (7  ) greater  than  unity.  The 
spectral  regions  of  breakdown  are  regions  of  small  absorption.  This  is  just 


J.  E.  A.  Selby,  Atr.no spheric  Transmittance  from  0.25  to  28. 5 pm: 
Computer  Code  LOWTRAN2,  AFCRL-72-0745,  Air  Force  Cambridge 
Research  laboratories,  Mass.  (29  December  1972). 
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Fig.  19.  Band  Model  Transmittance  Results  for  Slant  Path 
Using  LINAVE2  Parameters 
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Model  Transmittance  Results  for  Slant  Path  Using 
•al  Dynamics  Parameters 


the  condition  for  which  the  CG  approximation  was  predicted  (paragraph  II.  D) 
to  fail  most  noticeably. 

A comparison  of  the  band  model  and  line-by-line  results  for”  is 
given  in  Fig.  23  for  the  slant  path  and  Fig.  24  for  the  horizontal  path.  The 
agreement  is  not  quite  so  good  as  that  obtained  for  7 although  the  reasons 
fo:  the  disci  epancies  are  probably  the  same. 

The  results  described  so  far  pertain  to  a single  line  of  sight  that 
passes  perpendicularly  through  a diameter  position  of  the  equivalent  homo- 
geneous cylindrical  plume  Application  of  transmittance  results  to  the 
analysis  of  data  from  real  systems,  however,  would,  in  many  cases, 
require  the  average  results  for  the  effects  of  the  entire  plume.  For  the 
cylindrical  plume  being  considered,  there  is  no  axial  variation  of  the  plume 
properties  and  the  average  effect  of  the  whole  source  can  be  obtained  by  an 
analysis  of  the  effects  caused  by  moving  the  line  of  sight  off  the  diameter 
position.  Since  the  temperature  and  species  concentrations  are  also  homo- 
geneous, the  net  result  of  this  \ariation  is  to  shorten  the  length  of  the  hot 
source  path.  Results  of  such  an  off-axis  variation  at  v = 3700  cm”1  for  the 
slant  path  viewing  geometry  are  given  in  Fig.  25.  The  effect  of  displacing 
the  line  of  sight  off  axis  (decreasing  the  source  path  length)  is  to  make  7 
smaller;  the  value  of  at  the  diameter  position  is  U.  205  while  near  the 
edge  its  value  is  only  0.  100.  The  average  transmittance  of  the  atmosphere 
to  the  entire  source  is  0.  190,  which  is  a difference  of  only  about  6 percent 
from  the  value  for  the  diameter  position  alone.  The  transmittance  value  for 
the  atmosphere  alone  (7)  is  0.  456.  It  is  evident  from  this  result  that  the 
effects  of  line  correlation  between  emitter  and  absorber  become  more  impor- 
tant as  the  source  becomes  more  optically  thin.  This  tendency  is  consistent 
in  that,  as  the  optical  depth  of  the  source  is  increased,  line  broadening  and 
overlap  would  cause  the  source  radiation  to  become  more  and  more  continuum 
in  nature  and  t would  approach  t. 


(LS  Approximation) 


The  effect  of  line  correlation  on  the  transmittance  appropriate  to 
wide  spectral  bandpasses  was  determined  for  both  the  slant  and  horizontal 
paths.  ^Le  2.  7-jim  spectral  region  was  divided  (somewhat  arbitrarily)  into 
three  relatively  distinct  subregions.  The  bandpass  between  3350  and 
3575  cm  represents  the  low  wave  number  band  wing,  the  bandpass  between 
3575  and  3725  cm  represents  the  band  center,  and  the  bandpass  between 
3725  and  3950  cm  represents  the  high  wave  number  band  wing.  In  each 
of  these  bandpasses,  the  average  values  of  7 and  7g  were  computed  from 


T = 


t (v)  dv 


and 


T 

e 


R (v)  dv 

sensor 


// 


R 

source 


(v»)  dv> 


respectively.  The  integration  limits  and  i-  are  the  lower  and  upper  wave 

number  values  for  the  bandpass  respectively,  and  R and  R are 

sensor  source 

the  sensor  and  source  radiance  spectra  respectively.  For  continuum 

emitters,  the  expression  for  T reduces  to  that  for  T if  R does  not 

e source 

vary  much  with  v between  v l and  v^.  The  results  for  T and  T along  with 
the  ratio  T/T^  are  listed  in  Table  8.  For  both  the  slant  and  horizontal  paths, 
the  largest  discrepancy  between  T and  T ^ occurs  in  the  band  center  region 
although  the  effect  is  nearly  as  great  in  the  high  wave  number  region.  In 
the  latter  region,  the  emission/transmission  phenomena  are  controlled 
almost  exclusively  by  the  FLO  content  of  the  optical  path,  and,  although  the 
absolute  amount  of  water  vapor  in  either  of  the  paths  is  small,  the  effect 
of  line  correlation  is  pronounced  (e.  g.  , a factor  of  nearly  3 for  the 
horizontal  path). 
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c.  INTEGRAL  CONVERGENCE  CRITERIA 

The  radiance  and  transmittance  calculations  for  atmospheric  paths 
presented  in  this  report  rely  on  a numerical  integration  over  both  the 
atmos  pheric  portion  and  the  source  portion  of  the  total  optical  path.  The 
principal  integrations  required  are  those  necessary  to  compute  the  path- 
averaged  variables  x£  and  ^ [ Eqs.  (16)  and  (17)  respectively]  , the  trans- 
mittance t.  in  the  LS  approximation  [ Fq.  (26)]  , and  the  radiance  L[  Eq.  (7)] 
for  both  the  source  path  and  the  total  path.  All  of  these  integrations  are 
performed  using  a modified  Simpson  integration  routine.  The  principal 
variable  contributing  to  the  accuracy  of  the  integrations  is  the  number  N 
of  intervals  into  which  the  integration  oath  is  divided.  Present  calculations 
have  been  executed  by  dividing  the  integration  paths  into  equal  size  geometric 
intervals.  The  general  criterion  for  computational  accuracy  is  that  the  size 
of  the  integration  interval  As  should  be  smaller  than  or  of  about  the  same 
magnitude  as  the  scale  length  for  significant  change  in  the  variables  being 
integrated. 

The  calculations  for  and  p^  require  evaluation  of  the  two  integrals 

/s  _ 

k(s')  c(s')  p(s')  ds' 
u 

I?(b)  = f k(s')  c(s')  p(s')  p(s')  ds' 

J0 

If  the  atmosphere  is  approximated  as  plane-parallel,  isothermal, 
exponential  and  if  p is  written  as  P = pQp,  then 

I.(s)  = kp.  rS  c(s‘)e-s'cos  e/hdS' 

■'o 
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I2(s)  = k P0Pl  X C(S,)  e 


-28'  cos  G/h 


where  Pj  is  the  pressure  at  s = 0,  G is  the  zenith  angle  to  the  positive  s 
direction,  and  h is  the  exponential  scale  height.  For  the  case  where  c(s)  is 
nearly  constant  (e.  g.  , for  C02),  the  scale  lengths  for  significant  variations 
of  the  integrands  of  Ij  and  I2  are  h secG  and  (h/2)  secG  respectively.  In 
order  to  satisfy  the  more  stringent  of  the  two  convergence  criteria  implied, 
the  condition  As  ~ (h/2)  secG  should  be  met.  The  path  length  in  a plane 
parallel  atmosphere  between  altitudes  z^  and  z2  (z2  > z^)  is  simply 
S = (z  - z ) secG.  The  numl  er  of  intervals  into  which  the  path  should  be 
divided  to  satisfy  As  £ (h/2)  secG  is  thus 


2(z,  - z.) 


For  most  of  the  calculations  presented  here,  the  atmosphere  is  assumed  to 
end  at  an  altitude  of  100  km.  For  the  75-deg  zenith  angle  path  oi  iginating  at 
z = 20  km,  then,  N should  be  greater  than  about  25  if  a scale  height  of 
h cr  6.  5 km  is  assumed. 

For  the  case  v/here  c(s)  is  not  constant,  a scale  length  can  still  be 

estimated  by  approximating  c(s)  in  an  exponential  form.  For  example,  in 

7 

the  first  12  km  of  the  tropical  model  atmosphere  of  McClatchey  et  al.  , 


c(s)  for  H2<3  can  be  approximated  as 


cH2O(z|=cH2O(0)e 


-z/2.  17  km 


Thus,  the  scale  length  for  the  integrand  of  I2  becomes  (2/h  + 1/2.  17  km) 
secG.  For  slant  paths  originating  near  ground  level  (z^  0),  the  scale 


length  then  is  about  1.  7 km  and  N should  be  greater  than  about  77.  More 
rapid  variations  in  cH^Q  are  sometimes  present  at  higher  altitudes,  but 
tile  absolute  value  of  is  so  low  that  little  error  can  result  by  not 

integrating  through  these  regions  with  the  required  interval  size.  For 
horizontal  (6  = 90  deg)  atmospheric  paths,  the  scale  lengths  become  infinite 
and  imply  that  no  integration  is  necessary. 

Integration  of  Eqs.  (26)  and  (7)  involves  spatial  variation  of  the 
transmittance  derivative.  The  local  scale  length  for  variation  of  this 
parameter  is  (k  c p)  . For  atmospheric  paths,  k can  change  greatly  with 
wave  number,  and  c and  p can  again  change  with  path  position.  A worst 
case  estimate  of  the  required  number  of  intervals  required  for  an  accurate 
integration  can  be  obtained  by  computing  the  largest  value  ofkcp  occurring 
along  a path.  For  slant  paths,  this  occurs  at  the  lower  altitude  level  where 
p is  largest  (anr  usually  c also)  and  for  the  spectral  position  of  strongest 
absorption.  l or  example,  for  the  75-deg  slant  path  originating  at  20  km, 
the  conditions  for  CO^  absorption  are  paO.  0558  atm,  c^.^  cr  0.  00033, 
k (207  K,  3700  cm  ) ~ 1.10  cm  */atm  so  that  (k  c p)  * 0^49  km.  Thus, 

-he  number  of  intervals  should  he  about  160.  For  the  100-km-long  hori- 
zontal path  situated  at  20-km  altitude,  N should  be  greater  than  200.  The 
value  of  N = 160  for  the  slant  path  is  a worst  case  estimate,  while  the  value 
of  N = 200  for  the  horizontal  path  is  near  accurate. 

The  actual  convergences  displayed  by  numerical  calculations  using 
ATLES  are  presented  in  Figs.  26  and  27  for  the  slant  and  horizontal  paths 
considered  thus  far.  For  these  calculations,  values  of  N for  the  atmospheric 
path  (NAPATH)  of  10,  20,  30,  50,  100,  150,  and  200  (the  limit  of  the  com- 
puter program)  were  used.  For  the  source  portion  of  the  path,  N was  held 
constant  at  NSPATH  = 30.  This  value  of  NSPATH  is  easily  sufficient  to 
provide  less  than  1 percent  accuracy  for  calculations  within  the  source. 

For  the  slant  path  (Fig.  26),  both  t and  are  seen  to  converge  at 
NAPATH  30.  This  value  is  about  what  would  be  expected  on  the  basic  of 
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the  criterion  required  for  the  calculation  of  xg  and  (—25),  but  is  much 
less  than  that  estimated  for  the  convergence  of  the  radiance  and  trans- 
mittance equ^ions  (-160).  This  result  confirms  the  worst  case  nature  of 
the  estimate.  For  the  horizontal  path,  a quite  different  convergence 
characteristic  is  evident  (Fig.  27).  The  results  for  7 and  ~e  using  the  CG 
approximation  are  accurate  for  even  the  smallest  NAPATH  value  used 
(NAPATH  = 5)  while  the  results,  at  least  for  t,  using  the  LS  approximation 
are  still  in  error  by  more  than  20  percent  at  the  largest  value  of  NAPATH 
used  (NAPATH  = 200).  These  results  are  explainable  with  the  estimations 
presented  earlier.  For  a horizontal  path,  little  consideration  has  to  be 
given  to  the  convergence  of  and  pe  since  the  path  is  nearly  homogeneous 
and  th evaluation  of  these  parameters  simply  involves  multiplication.  Also, 
in  the  CG  approximation,  the  transmittance  is  given  by  an  analytical  function 
( Eq.  (l8)]so  that  no  integration  is  involved  here  either.  Only  for  the 
numerical  calculation  of  the  transmittance  in  the  LS  approximation  | Eq.  (26)] 
is  an  integration  involved.  The  previous  estimate  for  this  integration  to  be 
accurate  at  3700  cm-1  was  in  fact  NAPATH  > 200.  Thus,  the  estimation 
technique  is  confirmed. 

A serious  problem  is  presented  here.  For  calculations  over  very 
long,  nearly  horizontal  atmospheric  paths  within  the  LS  approximation,  the 
number  of  intervals  into  which  the  path  should  be  divided  to  obtain  aicui  ~te 
results  may  be  prohibitively  large.  The  solution  used  in  the  present  work 
has  been  to  dispense  with  the  calculation  of  the  strictly  atmospheric  trans- 
mittance properties  in  the  LS  approximation  and  use  only  the  results  of  the 
CG  approximation.  This  can  be  done  with  quite  good  accuracy  since  the 
degree  of  inhomogeneity  of  the  atmosphere  is  not  too  large  to  be  handled 
adequately  by  the  CG  approximation.  Results  obtained  for  various  atmos- 
pheric paths  in  both  the  CG  and  LS  approximations  (and  under  conditions 
where  the  results  of  the  latter  approximation  were  confirmed  to  be  accurate) 
have  never  shown  a difference  of  more  than  1 percent  between  the  two 
approximations.  Calculations  using  the  LS  approximation  hav*  been  excluded 
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when  necessary  in  the  computation  of  strictly  atmospheric  radiance/ 
transmittance  results,  but  the  approximation  is  always  retained  for  compu- 
tations involving  the  whole  (atmosphere  plus  hot  source)  optical  path. 


VI.  SUMMARY 


A.  SUMMARY  OF  PRESENT  WORK 

A band  m jdel  formulation  designed  for  computing  the  effective 
transmittance  uf  the  atmosphere  to  hot  H^O/CO^  plume  emission  sources 
has  been  prese  ited  and  discussed  in  detail.  This  model  is  formulated 
within  the  con"  nes  of  the  statistical  band  model  for  a random  array  of 
pressure-bri.  i ' ined  Lorentzian  lines  of  equal  strength  and  width.  The 
effects  of  line  correlation  between  the  plume  emission  spectrum  and  the 
atmospheric  H^O/CO,  absorption  spectrum  are  accounted  for  by  treating 
the  entire  optical  path  extending  from  the  sensor  position  through  the 
intervening  atmosphere  and  including  the  hot  emission  source  as  a single 
radiating/ absorbing  entity.  The  high  degree  of  inhomogerr ity  and  non- 
isothermality  of  such  an  optical  path  are  accounted  for  by  either  the  Curtis- 
Godson  or  Lindquist-Simmons  approximation.  In  addition,  the  inhomo- 
geneity conditions  for  which  the  Curtis- Godson  approximation  is  likely  to 
fail  are  determined.  Procedures  are  presented  for  the  efficient  application 
of  the  Lindquist-Simmons  approximation  (which  does  not  display  the  faults 
of  the  Curtis-Godson  approximation). 

Calculations  within  this  formulation  and  appropriate  to  the  viewing 
geometries  of  satellite  surveillance  systems  and  aircraft-to-aircraft  field 
measurements  programs  have  been  performed.  For  selected  wide  spectral 
bandpass  regions,  the  effect  of  line  correlation  is  shown  to  result  in  dis- 
crepancies of  as  much  as  a factor  of  3 between  currently  calculated  and 
correctly  calculated  average  transmittance  values. 

B.  FUTURE  WORK  ARE  *S 

The  course  of  the  present  work  has  indicated  three  general  areas  in 
which  future  work  would  be  highly  desirable.  The  first  area  concerns  the 
band  model  formulation  itself.  This  model  was  developed  within  the 
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statistical  band  model  of  equal- strength  Lorentzian  lines  primarily  so  that 
the  Lindquist-Simmons  approximation,  as  developed  by  Lindquist  and 
Simmons  , could  be  incorporated  directly.  Although  the  statistical  model 
is  applicable  to  the  consideration  of  absorption  by  H^O,  its  application  to 
C02  may  not  be  as  accurate  as  the  application  of  a regular-line  band  model. 
Work  on  the  development  of  an  approximation  analogous  to  the  Lindquist- 
Simmons  approximation  but  applicable  to  band  models  other  than  the  statis- 
tical model  for  equal -intensity  Lorentzian  lines  is  thus  indicated.  A 
generalization  to  line  shapes  other  than  simple  Lorentzian  would  also  be 
highly  desirable. 

The  second  area  concerns  the  band  model  parameters  k and  PQ  used 
with  the  band  model.  Adequate  band  model  parameter  sets  exist  for  the 
high- temperature  regions  applicable  to  combustion  gas  plumes  and  for  the 
low-temperature  regions  applicable  to  atmospheric  absorption  paths.  If  the 
problem  of  atmospheric  transmittance  could  be  treated  in  a manner  in  which 
the  source  emission  and  atmospheric  absorption  were  discoupled  processes, 
these  existing  parameter  sets  would  all  that  is  needed  to  handle  the 
transmittance  calculations.  Correct  computation,  however,  requires  the 
entire  optical  path  to  be  treated  as  a coupled  system.  Thus,  a consistent 
set  of  band  model  parameters  applicable  to  the  entire  temperature  range 
from  atmospheric  temperatures  to  combustion  gas  temperatures  is  required. 
Here,  the  term  "consistent"  implies  not  only  that  a realistic  temperature 
variation  for  k and  Pq  be  found,  but  also  that  these  parameters  have  the 
same  meaning  at  all  temperatures.  For  example,  k and  p^  derived  from 
the  AFCRL  line  compilation  have  explicit  meanings  by  virtue  of  theii  specific 
defining  equations  (Eqs.  [43]  and  [ 50 ] respectively).  They  are,  respectively, 
the  mean  absorption  coefficient  in  Av  and  the  mean  line  width  to  spacing 
ratio  in  Av  (for  nonresonant  self-broadening).  The  quantities  from  the 
General  Dynamics  tabulation,  at  least  for  H^O,  which  bear  the  same 
symbolic  names  are,  however,  simply  empirically  measured  constants 
which  best  fit  the  band  model  representation  chosen  to  analyze  the  data. 
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although  the  correspondence  of  the  General  Dynamics  k to  the  mean 
absorption  coefficient  and  the  General  Dynamics  Pq  to  the  mean-line-width- 
to- spacing  ratio  is  probably  quite  close.  The  fact  remains  that  a set  of 
experimentally  derived  k and  p^  values  can  yield  highly  accurate  trans- 
mittance results  and  yet  not  correspond  to  their  intu’tive  physical  meaning. 
Hence,  the  generation  of  consistent  band  model  parameter  sets  for  H-,0 
and  CO ^ that  cover  the  entire  temperature  range  required  for  the  calcula- 
tions of  atmospheric  transmittance s to  hot  gas  emission  sources  should  be 
given  serious  consideration. 

The  third  general  area  is  simply  one  of  completeness.  The  present 
work  has  considered  only  the  emission /absorption  effects  of  CC^  and  H^O 
within  both  the  source  and  the  absorbing  atmosphere.  The  effects  of 
radiation  and  absorption  by  minor  molecular  constituents  as  well  as  the 
effects  of  particulate  emission  in  the  source  (e.g.,  carbon)  and  particulate 
extinction  (absorption  or  scattering)  in  the  atmosphere  should  be  included 
to  give  a more  complete  description  of  the  problem.  Consideration  of 
background  signals  (e.g.,  earth,  sky,  and  missile  engine  hot  parts)  should 
also  be  undertaken 
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LABORATORY  OPERA  1 IONS 


The  Laboratory  Operations  of  The  Aeroapace  Corporation  ia  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encounteri  j in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
senslti-  a materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation -t  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Research  Laboratory:  Electromagnetic  theory,  dei'ces,  and 
propagation  phenomena.  Including  plasma  electromagnetics,  quantum  electronics, 
lasers,  and  electro-optics:  communication  sciences,  applied  electronics,  semi-’ 
conducting,  superconducting,  and  crystal  device  physics,  optical  and  acoustical 
Imaging  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon,  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  anti  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Physics  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  ol  the  atmosphere,  aurorae 
and  alrglow,  magnetospheric  physics,  uosmic  rays,  generation  and  propagation 
of  plasma  waves  In  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields,  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth’s  atmosphere,  ionosphere,  and 
magnetosphere,  the  effects  of  optical,  electromagnetic,  and  particuiate  radia- 
tions in  space  on  space  systems. 
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